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The  topic  of  this  report  is  cavita_iqon  nucleation  and 

concerns  itself  With  both  the  cavitation  nucleus  and  the  mechan 
isms  by  which  the  nucleus  is  stabilized.  /Three  different. theor 
ies  of  ttie  cavitation  nucleus  are  .eunsidered  most  plausible. 

They  are  'the  crevice  model ,  the  surfactant  skin  model,  and  the 
ionic  skin  model.  ^  These  hypothesized  nuclei  are  physically 
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different  and  each  has  been  used  to  explain  separate  groups 
of  data.  Presented  in  this  report  are  theoretical  and  experi¬ 
mental  advances  made  while  searching  for  a  unified  model  of  a 
cavitation  nucleus. 


.The  previous  crevice  model  ‘Trail''  r  eq  u  i  r  e  d  thatj  in  order  for 
a  vaporous  cavity  to  grow  from  a  crevice  nucleus,  the  liquid-  . 
gas  interface  must  reach  the  receding  contact  angle.  This 
condition  is  necessary/but  not  sufficient.  An  additional 
criterion  is  that  the  radius  of  curvature  of  the  interface  must 
be  greater  than  a  critical  radius,  beyond  which  the  cavity  is 
mechanically  unstable.  Using  these  nucleation  criteria,  the 
model  is  rederived. ^3  The  results  of  the  rederivation  are  quite 
satisfactory.  ^Both  -the>old  and  new  crevice  models  predict  the 
cavitation  threshold  as  a  function  of  surface  tension  and  temp¬ 
erature  well  .  However,  they feo dolsldiverge  when  predictions  of 
the  threshold  as  a  function  of  dissolved  gas  content  are  consid 
ered.  The  old  cr#vice  model  predicts  a  linear  dependence  but 
experiment  shows  -tha-t,  the  threshold  increases  rapidly  with 
decreasing  gas  content.  This  behavior  is  predicted  by  the 
revised  model. 

,  •  "  ’  ;  -  c  ■  V a.  — 

In  addition  to  vaporous  cavitation,  the  revised  model  is 
applied  to  diffusion  cavitation^  ^--ia  shown  that  the  model'--, 
-predicts  features  present  in  the  data  previously  explained  only 
in  terms  of  the  var  y  ing-permeab  il  ity  model  . 

The  experimental  part  of  the  dissertation  consists  of  a 
series  of  experiments  performed  to  measure  the  acoustic  cavita¬ 
tion  threshold  of  water  as  a  function  of  the  dissolved  ion 
concentration.  The  results  indicate  that  up  to  about  1  mrool/1, 
the  threshold  increases  with  increasing  concentration.  This 
is  directly  opposite  to  predictions  based  on  the  ionic  skin 
model.  However,  the  results  of  these  experiments  are  difficult 
to  explain  with  any  of  the  models  in  their  present  state  of 
development . 
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Introduction 


1-1  Statement  of  the  Problem 

The  topic  of  this  dissertation  is  stabilization  mechanisms  for  the 
cavitation  nucleus.  For  the  purposes  of  this  dissertation,  cavitation 
may  be  defined  as  the  formation  of  vapor  cavities  or  gas  bubbles  in  a 
liquid  through  a  reduction  of  the  pressure  within  that  liquid.  Although 
the  majority  of  this  work  will  pertain  to  acoustic  cavitation,  many 
types  of  cavitation  are  not  acoustic  in  origin.  These  processes  will 
not  be  ignored.  Cavitation  occurs  in  a  wide  variety  of  placet,  ranging 
from  sterile  laboratory  systems  to  dirty  beer  mugs,  from  ship  propellers 
to  inside  biological  systems.  That  cavitation  occurs  in  these  vastly 
different  environments  can  be  ascribed  to  a  single  common  phenomenon:  a 
(microscopic)  site  exists  at  a  place  within  the  system  where  conditions 
are  favorable  for  the  formation  of  a  (macroscopic)  bubble.  This  site  is 
called  a  cavitation  nucleus. 

The  impetus  behind  this  research  was  to  answer  the  questions:  what 
is  the  nature  of  this  nucleus?  Is  a  nucleus  unique  to  a  given 
environment,  or  can  a  nucleus  cause  cavitation  in  more  than  one 
environment  and  under  more  than  one  set  of  circumstances?  indeed,  can 
all  cavitation  processes  be  attributed  to  a  single  type  of  nucleus?  The 
answers  to  these  questions  have  been  sought  for  more  than  a  century  and 
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completely  satisfactory  answers  will  not  be  given  here.  In  a  practical 
sense,  it  was  suspected  at  the  outset  of  this  project  that  all  of  the 
answers  would  not  be  found.  They  were  meant  to  serve  as  "the  carrot  in 
front  of  the  horse".  What  will  be  presented  here  are  some  theoretical 
and  experimental  advances  made  while  seeking  the  answers.  Not  only  do 
these  results  clarify  some  of  the  past  and  present  research  in  this 
field,  they  also  serve  as  beacons  pointing  to  new  directions  that  future 
researchers  should  pursue. 

1—2  Historical  Perspective 

Perhaps  the  first  type  of  cavitation  observed  was  the  formation  of 
bubbles  in  liquids  supersaturated  with  gas.  This  is  familiar  to  all  of 
those  who  have  ever  opened  or  poured  a  carbonated  beverage.  Tomlinson5 
discussed  a  series  of  twelve  experiments  performed  with  soda-water  and 
various  solids.  He  contended  that  if  the  solids  were  "chemically 
clean",  no  bubbles  would  form  on  them.  He  went  on  to  point  out  that  If 
a  solid,  initially  chemically  clean,  came  into  contact  with  a  dirty 
cloth,  unclean  air,  or  dust  particles,  copious  bubbles  formed  on  it  when 
immersed  in  the  soda-water.  Tomlinson  concluded  that  chemically  clean 
solids  are  perfectly  wetted,  whereas  dirty  solids  are  not.  From  this  he 
argued  that  if  a  solid  were  dirty  it  would  have  a  s .  ■■  \  r  adhesion 
force  for  the  gas  then  for  the  liquid  and  that  it  wa  .  asymmetry 
that  caused  bubble  formation.  Although  this  line  of  reasoning  is 
incorrecc,  the  results  of  his  extensive  experimentation  with 
effervescence  support  the  explai nation  offered  by  his  contemporary  M. 
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Gernez  .  Gernez  proposed  that  the  outgassing  of  supersaturated  gaseous 
solutions  was  caused  by  gas  pockets,  embedded  in  the  fissures  found  on 
the  surface  of  every  solid,  regardless  of  its  degree  of  finish.  (It  is 
therefore  Gernez  who  is  the  initial  proponent  of  the  crevice  model  of 
the  cavitation  nucleus  -  a  model  to  be  discussed  extensively  here.  ) 
Effervescence  is  a  form  of  cavitation  which  requires  gaseous  diffusion 
to  "pump  up"  gas-containing  nuclei  untii  they  grow  large  enough  to  rise 
to  the  surface  of  the  liquid.  At  about  the  same  time  that  Tomlinson  and 
Gernez  were  doing  thair  work,  another  genre  of  cavitation  was  being 
investigated. 

This  other  kind  of  cavitation  has  its  origin  with  the  conception 

3 

and  development  ot  hydrodynamics.  Not  until  Daniel  Bernoulli  set  down 

the  guidelines  for  this  new  branch  of  science  in  1738  was  it  understood 

that  a  negative  pressure  could  be  produced  in  a  liquid.  In  the 

following  years  Euler  and  d'Alembert  debated  over  the  consequences  of 
4 

negative  pressures  .  Euler  believed  that  the  application  of  a  negative 
pressure  could  result  in  a  rupture  of  the  liquid  while  d'Alembert 
refused  to  accept  this  view. 

The  rise  of  cavitation  as  a  topic  for  scientific  research  began 
with  the  development  of  high-powered,  high  rpm  steam  turbines  in  the  raid 
1800's.  With  this  carae  the  means  of  moving  an  object  (such  as  a 
propeller)  through  a  fluid  rapidly  enough  so  that  the  object  lost 
contact  with  the  fluid.  The  most  obvious  concern  to  which  this  new 
phenomenon  drew  attention  was  that  the  state-of-the-art  propeller  design 
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was  poor. 


A  second,  and  as  was  soon  evident,  much  greater  concern  was 


propeller  erosion.  In  fact,  this  problem  was  so  critical  that  in  1915 
the  British  Admiralty  appointed  a  special  subcommittee  to  investigate 
this  phenomenon.  In  1917,  Lord  Rayleigh"1  became  involved.  His  solution 
of  the  equations  governing  vapor  cavity  dynamics  is  still  quite  useful 
today. 

Concurrent  to  this  work  the  first  attempts  to  measure  the  tensile 
strengths  of  liquids  were  being  undertaken.  In  about  1850,  Berthelot 
began  to  measure  the  tensile  strength  of  water  and  found  it  to  be  on  the 
order  of  50  bars  (maximum)^.  His  method  consisted  of  heating  a  liquid, 
which  was  sealed  in  an  evacuated  glass  capillary  tube,  until  the 
expanding  liquid  completely  filled  the  tube.  The  temperature  at  which 
the  cube  was  filled  was  recorded  and  the  liquid  was  then  aLlowed  to 
cool.  Because  the  liquid  adhered  to  the  sides  of  the  tube,  this 
subsequent  contraction  resulted  in  the  liquid  undergoing  a  tensile 
stress,  eventually  ending  in  the  rupture  of  the  liquid.  The 
temperature  at  which  rupture  occurred  was  also  recorded.  From  a 
knowledge  of  the  filling  and  rupturing  temperatures,  along  with  the 
thermal  expansion  coefficient  for  the  liquid  being  tested,  Berthelot  was 
able  co  calculate  a  tensile  strength.  Many  other  measurements  of  the 
tensile  strength  of  water  through  these  static  means  have  been  made  over 
the  intervening  one  hundred  and  thirty-four  years,  and  their  results 
(ranging  from  as  low  as  about  5  Co  as  high  as  about  200  bars)  are 
comparable  Co  Berthelot's.  A  detailed  discussion  of  these  works  will 
not  be  given  here.  The  point  is  that,  for  the  first  time  investigators 
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were  Intentionally  putting  liquids  under  stress  and  considering  the 
consequences.  The  application  of  a  tensile  stress  to  a  liquid  is  the 
basis  for  present  day  acoustic  cavitation  studies,  which  are  actually 
studies  of  the  "dynamic"  tensile  strength,  as  opposed  to  these  earlier 


"static"  measurements. 


Acoustic  cavitation  was  first  observed  during  the  period  from  IS  15 
to  1920  by  Langevin  and  his  co-workers  while  he  was  pioneering  the  field 
of  ultrasonics^.  Sb'llner®  was  one  of  the  first  to  observe  cavitation  in 
degassed  liquids  at  room  temperature  and  atmospheric  pressures  (1936). 
In  the  1930's  and  1940's  a  number  of  other  researchers  investigated 
various  aspects  of  ultrasonic  cavitation,  among  them  R.W.  Boyce** ,  E.N. 
Harvey*’,  and  Kc~nfcld  and  Suvorov®. 

The  fundamental  problem  of  cavitation  was  by  now  well  formed. 
Numerous  measurements  of  the  cavitation  threshold  (tensile  strength)  of 
water  yielded  results  in  the  range  of  5-25  bars  (using  the  most  careful 
and  painstaking  procedures  thresholds  of  up  to  300  bars  have  been 
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obtained  ,  but  they  are  unusual)  and  yet,  theoretical  predictions  of  the 
homogeneous  (pure  liquid)  threshold  of  water  are  thousands  of  bars.** 
What  accounts  for  this  order-of-magnitude  (best  case)  discrepancy? 
Failure  in  solids  under  teusion  is  usually  attributed  to  an  imperfection 
in  the  solid.  The  same  logic  (i.e.,  imperfections  in  the  liquid)  was 
applied  to  the  failure  of  liquids. 

In  the  remainder  of  this  section  different  types  of  cavitation  will 
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of  ultrasonics7.  Sb'llner^  was  one  of  the  first  to  observe  cavitation  in 

degassed  liquids  at  room  temperature  and  atmospheric  pressures  (1936). 

In  the  1930 '8  and  1940's  a  number  of  other  researchers  investigated 

various  aspects  of  ultrasonic  cavitation,  among  them  R.W.  Bcyce6,  E.N. 

Harvey  ,  and  Kornfeld  and  Suvorov  . 

The  fundamental  problem  of  cavitation  was  by  now  well  formed. 

Numerous  measurements  of  the  cavitation  threshold  (tensile  strength)  cf 

water  yielded  results  in  the  range  of  5-25  bars  (using  the  most  careful 

and  painstaking  procedures  thresholds  of  up  to  300  bars  have  been 

g 

obtained  ,  but  they  are  unusual)  and  yet,  theoretical  predictions  of  the 
homogeneous  (pure  liquid)  threshold  of  water  are  thousands  of  bars.^ 
What  accounts  for  this  order-of-magnitude  (best  case)  discrepancy? 
Failure  in  solids  under  tension  is  usually  attributed  to  an  imperfection 
in  the  solid.  The  same  logic  (i.e.,  imperfections  in  the  liquid)  was 
applied  to  the  failure  of  liquids. 

In  the  remainder  of  this  section  different  types  of  cavitation  will 
be  mentioned  in  an  overview  of  the  field  in  the  last  half-century.  The 
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investigators*  of  these  branches  of  cavitation  research  sought  these 

imperfections  (cavitation  nuclei)  to  explain  the  various  aspects  of  the 

particular  cavitation  process.  Cavitation  nucleated  by  these 

inhomogeneities  is  termed  heterogeneous  cavitation. 

In  1944,  Dean^  suggested  that  vapor  cavities  form  in  the  center  of 

small  vortices  produced  by  turbulent  motion  of  water  around  solid 

objects.  In  the  same  year,  Harvey^  suggested  that  gas  pockets 

stabilized  at  the  bottom  of  crevices  found  on  dirt  particles  accounted 

for  the  cavitation  found  ir>  animals.  He  had  concerned  himself  mainly 

with  cavitation  formed  by  c.ompression/decotnpression  processes  in 

animals  (the  "bends");  however,  he  did  discuss  his  model's  application 

to  acoustic  cavitation.  One  of  the  most  careful  and  extensive 

12 

experiments  of  the  time  was  performed  by  Briggs,  Johnson,  and  Mason 
They  measured  cavitation  thresholds  for  a  variety  of  liquids  as  a 
function  of  dissolved  gas  pressure,  viscosity,  and  pulse  length. 
However,  their  theoretical  explanation  of  their  results  was  later  shown 
to  be  incorrect^. 

The  1950's  was  a  decade  of  diversification  for  cavitation  research. 

Perhaps  the  most  important  example  of  this  breadth  is  the  invention  of 

the  bubble  chamber  by  Glaser.  Tnis  invention  gave  rise  to  the  field  of 

1 3 

radiation-induced  cavitation.  Seitz  gave  an  explanation  for  bubble 
formation  in  superheated  liquids  by  charged  particles.  He  concluded 
that  the  majority  of  bubbles  were  "nucleated  by  moderately-energetic 
free  electrons  produced  by  the  incident  particles  in  Coulomb 
encounters".  The  bubble  chamber  was  developed  as  a  tool  for  elementary 
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particle  phvsics  while  radiation-induced  cavitation  research  took  a 

1 A 

different  avenue.  Lioberman  measured  the  threshold  for  acoustic 

cavitation  in  pentane  and  acetone  exposed  to  neutron  and  beta  sources. 

He  concluded  that  the  cavitation  was  nucleated  by  the  recoil  of  carbon 

ions;  these  results  agreed  well  with  Seitz's  theory. 

In  the  field  of  acoustically  induced  cavitation,  a  different  type 

1 5 

of  nucleus  was  suggested  in  1954  by  Fox  and  Herzfeld  .  This  nucleus 

consisted  of  a  gas  bubble  surrounded  by  a  rigid  skin  of  organic 

molecules.  This  organic  skin  model  will  be  discussed  further  in  the 

next  section.  An  extensive  set  of  measurements  was  reported  by 
1 6 

Strasberg  .  He  analyzed  his  results  in  terms  of  the  nuclei  proposed  up 


to  that  time  and  found  the  crevice  model  offered  the  best  explanation. 

Other  noteworthy  papers  of  the  195C's  are  by  Connolly  and  Fox^,  and 

Galloway^;  some  of  their  results  will  be  discussed  later. 

1  9 

Sette  and  Wanderlingh*  were  the  first  to  study  acoustic  cavitation 
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in  water  induced  by  cosmic  rays.  They  concluded  that  the  nuclei  were 

recoiling  oxygen  atoms  and  their  results  agreed  well  with  calculations 

20 

based  on  Seitz's  theory.  The  Soviet  researcher  Akulichev  proposed 
another  cavitation  nucleus  based  on  bis  measurements  of  the  influence  of 
dissolved  ions  on  cavitation  thresholds.  His  nucleus,  similar  to  that 
of  Fox  and  Herzfeld,  consisted  of  an  ionic  skin  surrounding  a  gas 
bubble.  Greenspan  and  Tschiegg^  described  a  method  bv  which  they  were 
able  to  attain  unusually  hieh  cavitation  thresholds.  This  method 


consisted  of  prolonged  filtering  and  cleaning  of  the  liquids  to  be 
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tested  in  a  closed  system. 

22 

In  1970,  Apfel  developed  the  crevice  model  further  to  include  the 

size  of  the  crevice;  he  suggested  that  there  existed  n  "critical"  size 

for  the  crevice.  Crevices  larger  than  and  smaller  than  this  critical 

size  exhibit  quite  different  behavior.  Apfei  published  papers  on  both 

23 

homogeneous  and  heterogeneous  cavitation.  In  1971  ,  he  modified  a 

theory  for  homogeneous  nucleation  of  vapor  cavities  at  the  interface 

between  a  flat  solid  and  a  pure  liquid  to  account  for  nucleation  at  the 

Interface  between  two  pure  liquids.  Apfel's  crevice  model  was  modified 
24 

further  by  Crum  who  considered  the  variation  of  contact  angles  with 

surface  tension.  The  crevice  model  in  this  latest  form  was  successful 

in  explaining  a  wide  variety  of  data  obtained  from  acoustic  cavitation 

25 

measurements.  WInterton  also  modified  Apfel's  model  and  used  it  to 

explain  boiling  and  non-acoustic  cavitation. 

At  about  the  same  time  Yount  was  developing  a  different  type  of 
26 

nucleus.  Modifying  the  organic  skin  models  of  Fox  and  Herzfeld,  and 
27 

Slrotyuk  ,  Yount  considered  his  nucleus  to  be  a  bubble  surrounded  by 

polar  surfactant  molecules.  Ihis  skin  has  the  property  of  offering  a 

varying  permeability  to  gas,  depending  upon  the  ambient  conditions.  He 

has  applied  his  model  to  corapression/decompresslon  processes  and  has 

been  successful  in  explaining  his  data. 

And  so  this  was  the  state  of  the  field  when  this  project  was 

28 

undertaken.  Since  that  time,  Yount  has  further  developed  his  model 
theoretically  and  experimentally  so  that  there  is  little  doubt  that  it 
is  plausible.  However,  the  crevice  model,  with  new  developments  to  be 


discussed  later,  explains  a  very  broad  set  of  data. 


After  this  brief 


and  sketchy  historical  review  of  the  field  in  general,  the  cavitation 
nucleation  and  stabilization  mechanisms  will  be  discussed  in  detail  in 
the  following  section  with  emphasis  on  the  organic  skin  model  of  Yount 
and  the  crevice  model  as  developed  by  Apfel  and  Crum. 

1-3  Review  of  Stabilization  Mechanisms 

The  most  obvious  candidate  for  the  cavitation  nucleus  is  a  gas- 

filled  bubble.  However,  this  nucleus  is  unstable.  A  bubble  will 

dissolve  in  an  undersaturated  solution  and  the  effect  of  surface  tension 

will  cause  it  to  dissolve  in  a  saturated  solution.  In  supersaturated 

solutions,  a  bubble  can  be  in  equilibrium  because  the  tendency  for  the 

bubble  to  dissolve  due  to  surface  tension  is  opposed  by  the  tendency  for 

the  bubble  to  grow  by  diffusion  of  gas  intc  it.  However,  this 

equilibrium  is  unstable;  the  bubble  will  grow  or  dissolve  depending  on 

whether  the  perturbation  increases  or  decreases  the  bubble's  radius 

29  30 

relative  to  its  equilibrium  radius.  ’  Therefore,  a  liquid  would  be 

free  of  bubbles  after  a  short  period  of  time  and  cavitation  thresholds 
would  quickly  approach  those  for  homogeneous  nucleation.  This  does  not 
imply  that  gas  bubbles  could  not  serve  as  cavitation  nuclei.  It  does 
imply,  however,  that  in  order  for  gas  bubbles  to  serve  as  cavitation 
nuclei,  they  must  be  stabilized  at  a  size  small  enough  to  prevent  their 
rising  to  the  surface  of  the  liquid,  yet  large  enough  so  that  they  will 
grow  when  exposed  to  negative  pressure  as  lew  as  a  few  bars.  In  other 


words,  a  stabilization  mechanism  must  exist  for  a  gas  bubble  before  it 


can  act  as  a  cavitation  nucleus. 

That  cavitation  nuclei  do  in  fact  contain  gas  has  been  borne  out 

in  numerous  experiments.  Some  measurements  show  that  the 

cavitation  threshold  is  strongly  dependent  on  dissolved  gas  pressure, 

2 1 

except  fo.  extremely  clean  water.  Also,  samples  which  undergo  a 
compression/  decompression  cycle  prior  to  testing  show  higher  thresholds 
than  samples  which  do  not. 

Various  types  of  stabilizing  skins  have  been  proposed.  These  skins 
usually  consist  of  contaminants  which  somehow  deposit  themselves  on  the 
bubble's  surface  and  counteract  the  surface  tension.  The  first  of  these 

skin  models  to  be  discussed  is  the  ionic  skin  model  proposed  by 

.  .  20 
Akulicnev. 


1-3.1  Ionic  Skin  Model 
31 

Alty  noted  that  it  has  been  known  since  the  1860's  that  liquid- 

gas  interfaces  are  electrically  charged.  He  experimentally  determined 

that  bubbles  quickly  acquired  a  charge  after  being  introduced  into  water 

and  behaved  as  though  they  carried  a  net  negative  charge.  In  a  later 
32 

paper"  ,  Alty  postulated  a  charging  mechanism  by  which  a  bubble  became 
charged  by  the  selective  adsorption  of  ions  onto  the  bubble  surface.  He 
noted  that  the  rate  at  which  the  bubble  dissolved  decreased  as  the 
charge  on  the  bubble  increased — a  phenomenon  which  he  attributed  to  the 
Coulomb  repulsion  of  the  adsorbed  charges.  The  mechanism  by  which  this 


selective  adsorption  took  place,  however,  was  unclear 
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Akulichev  postulated  that  not  only  could  the  rate  of  dissolution  of 
a  bubble  be  decreased,  the  dissolution  could  be  stopped  altogether.  He 
explained  the  selective  adsorption  in  terms  of  the  type  of  hydration 
exhibited  by  various  ions.  Those  exhibiting  positive  hydration  behave 
in  a  hydrophilic  manner.  They  are  surrounded  by  water  molecules  when 
dissolved  in  aqueous  solutions.  On  the  other  hand,  ions  exhibiting 
negative  hydration  behave  in  a  hydrophobic  manner  and  are  not  "as" 
surrounded  by  water  molecules  as  are  hydrophilic  ions.  Hydrophobic 
ions,  therefore,  are  more  likely  to  stay  on  the  surface  of  a  bubble  once 
they  encounter  it  than  are  hydrophilic  ions,  in  other  words — selective 
adsorption.  Hydrophilic  ions  include  Na+,  Li+;  examples  of  hydrophobic 
ions  are  Mg  ,  K  ,  F  ,  Cl  ,  and  I  . 

In  order  to  test  these  ideas,  Akulichev  measured  the  cavitation 
threshold  as  a  function  of  dissolved  ion  concentration.  Some  of  his 
results  are  reproduced  in  figure  1.  Figure  la  shows  the  variation  in 
threshold  with  concentration  for  three  salts:  LiOH,  Kl,  and  KBr-  Since 
OH  is  a  natural  constituent  of  water,  Akulichev  argued  that  it  could 
not  affect  the  bubble.  Lithium  is  hydrophilic.  Therefore,  the  presence 
of  LiOH  in  the  water  should  have  no  effect  on  the  threshold.  It  is  seen 
from  figure  la  that  the  threshold  is  independent  of  concentration  for 
LiOH.  The  ions  K  ,  I  ,  and  Br  are  all  hydrophobic  so  they  should  all 
be  adsorbed  by  a  bubble.  This  results  in  the  bubble  being  charged. 
Therefore,  because  of  Coulomb  repulsion  of  the  ions,  the  equilibrium 
radius  is  larger  than  for  an  uncharged  bubble.  This  adsorption  should 


Threshold  (bar)  Threshold  (bar) 


1  3/$  Gas  Content 

a 

LiOll 

K  I 

ih| 

-■  T"  ”  "O  - 

9; 

6.0* - * - - * - 

0  10  20  30 
Ion  Concentration  (mmol/1) 


Ion  Concentration  (mmol/1) 


Figure  1.  Threshold  versus  Dissolved  Ion  Concentration 
for  various  compounds  and  dissolved  gas  pressures.  Re¬ 
produced  from  Akulichev  (reference  20) . 
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produce  a  lower  threshold  for  aqueous  solutions  of  KBr  and  K1  than  of 
L10H.  Notice  that  the  effect  tends  to  saturate  at  a  concentration  of 
about  5  mmol/].  He  attributes  the  slightly  lower  threshold  (about  5% 
lower)  of  KBr  compared  to  that  of  K1  to  the  difference  in  mobilities  of 
the  two  ions  Br  and  I  .  Figure  lb  is  a  similar  graph  for  the  salts 
NaOH,  NaCl,  and  NaF.  In  this  instance,  the  positive  ion  (Na+)  is 
hydrophilic.  The  differences  in  the  three  curves  are,  therefore,  due  to 
the  different  natures  of  OH  ,  Cl  ,  and  F  .  Again,  the  threshold  is 
lower  for  solutions  of  the  hydrophobic  anions. 

Akulichev's  model  requires  that  os  the  ion  concentration  approaches 
zero,  the  thresholds  for  the  various  salt  solutions  approach  the 
threshold  of  LiOH  and  NaOH.  In  other  words,  the  threshold  must  increase 
with  decreasing  ion  concentration.  It  is  important  to  notice  two  things 
about  Akulichev's  data.  First,  the  lines  drawn  are  not  theoretical 
calculations.  They  are  best  fit  lines  to  his  experimental  data.  Second, 
no  measurements  are  made  for  ion  concentrations  below  about  1  mmol/1. 
The  thresholds  at  zero  concentration  are  extrapolated  from  thresholds 
measured  at  concentrations  greater  than  about  1  mmol/1.  We  have 
extended  Akulichev's  experimental  results  to  concentrations  below 
1  mmol/1.  These  measurements  wll]  be  discussed  in  detail  in  a  later 
chapter. 

Akulichev  and  Alty  considered  these  stabilizing  ions  to  be  trace 
contaminants  in  their  water.  It  is  known  that  water  of  all  but  the 
purest  of  samples  contains  such  contaminants.  Organic  molecules  are 
also  found  as  contaminants  in  water  and  they  have  been  employed  by 
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others  as  stabilizing  agents.  The  first  such  attempt  was  by  Fox  and 
Herz  fe  Id  ,  ^and  so  their  model  will  be  the  first  to  be  discussed. 


1-3.2  Rigid  Organic  Skin  Model 

Fox  and  Herzfeld'^  proposed  that  surface  active  organic  molecules 

could  form  a  rigid  skin  around  a  gas  bubble.  This  skin  would  be 

impermeable  to  gas  diffusion  and  would  be  mechanically  strong  enough  to 

withstand  moderate  hydrostatic  pressures.  A  nucleus  would  cavitate  when 

the  acoustic  pressure  reached  a  (negative)  level  at  which  the  skin  would 

tear  because  of  the  pressure  difference  between  the  gas  pressure  inside 

and  the  acoustic  pressure  outside.  According  to  this  model,  the  skin 

wourd  have  a  characteristic  crushing  pressure-  As  long  os  the 

hydrostatic  pressure  remained  below  this  crushing  pressure,  the  nucleus 

would  be  unaffected  by  Increases  in  hydrostatic  pressure.  If  the 

hydrostatic  pressure  exceeded  the  crushing  pressure,  the  skin  would 

crumble  and  the  nucleus  would  dissolve.  This  would  result  In  an 

increase  in  the  cavitation  threshold  since  the  nuclei  would  no  longer  be 
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present.  It  was  this  property  which  forced  Hei’zfeld  to  later  retract 
his  model,  Strasberg^  had  found  that  the  threshold  does  increase  with 
hydrostatic  pressure.  However,  the  variation  is  continuous;  there  is  no 
lower  limit  to  this  variation  as  predicted  by  the  rigid  skin  model. 


The  idea  of  an  organic  skin  has  been  developed  further  by 
2728  28 

others  ’  ;  the  most  viable  form  will  be  discussed  next. 
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1-3.3  Varying  Permeability  Model 

Rather  than  the  rigid  skin  of  organic  molecules  proposed  by  Fox  and 
Herzfeld ,  Yount26,28  has  developed  a  stabilization  theory  in  which  the 
dissolution  of  gas  bubbles  is  halted  by  a  non-rigid  organic  skin.  This 
skin  has  mechanical  strength  against  compression,  but  none  against 
tension.  13. is  postulated  nucleus  is  depicted  schematically  in  figure  2. 

This  nucleus  is  formed  in  the  following  manner.  A  gas  bubble  is 

introduced  into  a  liquid  (water,  say)  and  begins  to  dissolve.  While  the 

dissolution  takes  place,  surface  active  molecules  accuramulate  on  the 

surface.  Eventually,  within  minutes  or  seconds,  the  density  of  these 

surfactants  is  such  that  they  resist  the  collapse  of  the  bubble.  This 

is  presumably  a  Coulomb  interaction.  These  molecules  are  polar  (c. 

perhaps  even  carry  a  net  charge)  and  align  themselves  such  that  their 

polar  "heads"  face  outward  (toward  the  water)  and  their  tails  (typically 

hydrocarbons)  stick  inward.  Thus,  as  the  density  increases,  the 

separation  of  the  heads  becomes  small  enough  for  dipole  fields  to  become 

important.  Of  course,  one  would  expect  that  if  the  surfactants  were 

charged,  the  density  at  which  stabilization  occurs  would  be  lower 

-2  .  .  -3 

because  monopolar  fields  vary  as  r  rather  than  r 

Once  the  critical  density  is  reached,  the  radius  of  the  nucleus  may 
change  through  changes  in  the  number  of  molecules  on  the  skin.  The 
equilibrium  condition  for  this  nucleus  is  that  the  electrochemical 
potentials  of  the  skin  and  the  reservoir  must  be  equal.  The  reservoir 
is  a  thin  (perhaps  monomolecular )  layer  of  non-aligned  surfactants  which 
surrounds  the  skin  of  aligned  molecules.  This  equilibrium  condition  can 
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be  satisfied  by  the  accretion  or  deletion  of  molecules  from  the  skin  (or 

reservoir),  resulting  in  a  "large-scale"  change  in  radius.  Having 

adjusted  the  radius  in  this  manner,  subseauent  "small-scale"  changes  in 

radius  occur.  These  changes  involve  only  the  adjustment  of  the 

separation  of  adjacent  skin  molecules,  not  a  change  in  the  number  of 

molecules.  At  equilibrium,  the  nucleus  can  be  thought  of  as  two 

concentric  shells  of  negligible  thickness — the  outer  one,  the  reservoir, 

is  in  contact  with  the  liquid  and  the  inner  one,  the  skin,  is  in  contact 

with  the  gas.  So  far,  this  model  has  only  been  used  to  make  predictions 

about  compression/decompression  processes  during  which  gas  diffusion  can 

28 

plav  a  major  role.  However,  Yount  has  calculated  characteristic  times 

in  which  a  skin  molecule  can  be  accreted  or  deleted,  due  to  an 

inequality  of  electrochemical  potentials,  and  found  them  to  range  from 

10  to  10  seconds.  The  shortest  times  approach  the  period  of 

acoustic  sLgnals  used  in  ultrasonic  studies.  This  would  imply  that 

large-scale  changes  in  radius  may  be  able  to  keep  up  with  the  variations 

in  applied  pressure  associated  with  acoustic  cavitation. 

As  mentioned  above,  Yount  has  performed  only  long-time  scale 

experiments.  A  typical  experiment  will  now  be  described.  During  the 

experiment,  gelatin  samples  (surface  tension  around  50-55  dvn/cm)  are 

subjected  to  a  pressure  schedule  which  is  typically  like  the  one  shown 

in  figure  3.  The  sample,  about  0.A  ml  in  volume,  initially  at  a 

pressure  P  ,  is  rapidlv  compressed  to  a  pressure  P  .  "Rapidly"  means 
o'  m 

that  no  gas  diffuses  out  of  the  nucleus  during  the  compression.  The 
sample  is  held  at  P  long  enough  f  >r  the  dissolved  gas  pressure  in  the 
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sample  to  come  to  equilibrium  at  P^;  i.e.,  it  becomes  saturated.  The 

holding  time  is  typically  5.25  hours.  The  ambient  pressure  is  then 

reduced  rapidly  to  resulting  in  the  growth  of  some  of  the  nuclei. 

This  growth  occurs  through  gas  diffusion  from  the  now  supersaturated 

liquid  into  the  nuclei.  Several  minutes  are  allowed  for  the  nuclei  to 

grow  to  visible  size.  Then  the  total  number  of  visible  bubbles  are 

counted  and  correlated  with  the  pressure  schedule.  Yount  defines  the 

crushing  pressure,  pcrush*  as  pm-p0»  an^  ttie  supers£*-uration  pressure, 

P  ,  as  P  -P-.  One  manner  of  data  presentation  is  a  graph  of  P  vs. 
ss  mf  ‘  orss 

P  ,  for  lines  of  constant  number  of  bubbles.  Such  a  graph  is 
crush 

reproduced  in  figure  4. 

The  isopleths  (lines  of  constant  number  of  bubbles)  in  figure  4 

* 

show  a  gradual  decrease  in  slope  at  1pcrus{l=*^>  -p0’  segments  of  the 

lines  on  either  side  of  this  critical  point  are  more  or  less  linear. 

This  behavior  is  explained  as  follows.  For  modest  values  of  P  ,  the 

crush 

skin  remains  permeable  to  gas  diffusion.  This  enables  the  gas  inside 

the  nucleus  to  remain  in  equilibrium  with  the  dissolved  gas  in  the 

* 

gelatin.  However,  at  a  critical  ambient  pressure,  P  “^>0+pcrusp1  >  the 
skin  molecules  becosie  so  tightly  packed  that  the  skin  becomes 
essentially  impermeable  to  gas  diffusion.  This  prevents  the  interior 
and  exterior  of  the  bubble  from  remaining  in  equilibrium.  Therefore, 


the  gas  pressure  inside  the  bubble  will  be  less  than  the  gas  pressure  of 
the  saturated  gelatin,  resulting  in  a  lower  "effective"  Psg  during 
decompression.  Throughout  decompression  the  skin  is  always  permeable, 


"varying- 


indepeodent  of  P 


This  process  gives  rise  to  the  name 


permeablility". 

Perhaps  it  is  useful  to  point  out  that  for  acoustic  processes 

P  «P  and  that  the  holding  time  is  usually  too  short  for  any 
crush  ss  s  1 

significant  gas  diffusion  to  occur.  Other  aspects  of  this  elegant  model 
will  be  discussed  later. 

1-3.4  Crevice  Model 

The  crevice  model  is  fundamentally  different  from  the  other  models 

discussed.  Its  stabilization  mechanism  does  not  require  that  surface 

tension  be  nullified.  Rather,  it  uses  surface  tension,  combined  with 

geometrical  considerations,  to  stabilize  a  gas  pocket  at  the  bottom  of  a 

crevice.  It  is,  perhaps,  the  most  “worked  over"  of  the  models.  This  is 

partly  true  because  it  can  account  for  a  wide  range  of  experimental 

data,  whereas  the  other  models  can  account  for  only  limited  sets  of 

data,  most  often  that  of  the  model's  proponent.  Since  the  crevice  mciel 

has  undergone  several  revisions,  it  would  be  unreasonable  to  discuss  all 

of  them  in  detail.  Therefore,  the  model  will  be  discussed  briefly  in 

34 

the  form  reported  by  Crum  ;  this  model  will  then  be  revised  in  the  next 
chapter. 

The  crevice  model  assumes  that  gas  pockets  are  stabilized  at  the 
bottom  of  conical  cracks  or  crevices  found  on  hydrophobic  solid 
impurities  present  in  the  water.  The  essential  features  of  the  model 
are  depicted  in  figure  5.  When  the  liquid  is  saturated  with  gas,  th.v 
liquid-gas  interface  is  essentially  flat.  However,  when  the  liquid  is 
degassed,  the  interface  bows  toward  the  apex  of  the  crevice.  This 
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Figure  5.  Diagram  of  a  gas-filled  conical  crevice 
embedded  in  a  hydrophobic  solid  a)  in  the  degassed 
state  and  b)  in  the  nucleation  phase.  (After  Crum 
reference  34)  . 
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behavior  occurs  for  the  following  reason.  In  a  saturated  solution,  the 
dissolved  gas  pressure,  P^,  equals  the  pressure  which  the  liquid  exerts 
on  the  interface.  We  call  this  the  liquid  pressure,  P^,  and  define  it  to 
be  the  sum  of  all  the  pressures  present — hydrostatic  and  acoustic. 
Diffusion  maintains  the  gas  pressure  in  the  nucleus  at  the  dissolved 
gas  pressure.  Assuming  that  vapor  pressure  is  negligible  for  the  time 
being,  we  have  P  “P  for  saturated  solutions.  Since  there  is  no 

"  s 

pressure  difference  across  the  interface,  it  is  flat  and  the  Laplace 

pressure  (2Y  / r,  where  y  •*  surface  tension  and  r  *  radius  of  curvature 

of  the  interface)  is  zero.  When  the  liquid  is  degassed,  P  is  less  than 

8 

P^  and  the  interface  bows  toward  the  apex.  This  curvature  glres  rise  to 
a  Laplace  pressure  which  exactly  equals  the  difference  between  P^  and 
P^.  Since  all  pressures  are  now  balanced,  the  nucleus  does  not  dissolve 
and  remains  intact. 

Suppose  we  have  a  degassed  liquid  and  the  liquid  pressure  is 
increased.  In  response  to  this  condition,  the  interface  bows  inward 
more,  until  it  reaches  the  advancing  contact  angle.  At  this  point  any 
subsequent  taotion  of  the  interface  involves  the  entire  interface 
advancing  as  a  whole  toward  the  apex.  As  it  advances,  the  radius  of 
curvature  necessarily  becomes  smaller,  since  the  angle  of  the  interface 
measured  from  the  crevice  wall  is  now  fixed.  As  the  radius  of  curvature 
decreases,  the  Laplace  pressure  increases,  eventually  becoming  high 
enough  to  balance  the  increased  liquid  pressure. 

If  the  liquid  pressure  is  now  suddenly  decreased  becoming  negative. 


as  during  the  negative  portion  of  an  acoustic  cycle,  the  interface  bows 
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outward  and  may  reach  the  receding  contact  angle.  If  the  receding 

contact  angle  is  reached,  any  subsequent  motion  of  the  interface 

involves  the  entire  interface  receding  away  from  the  apex.  This 

recession  results  in  a  decrease  in  the  Laplace  pressure,  since  the 

radius  of  curvature  increases.  With  the  Laplace  pressure  decreased  the 

gas  pressure  inside  the  nucleus  is  retarded  even  less.  Previous 

16  22  24  25 

researchers  *  ’  ’  assumed,  therefore,  that  the  criterion  for 

nucleation  from  a  crevice  was  that  the  interface  must  reach  the  receding 
contact  angle.  They  argued  that  once  the  interface  had  reached  the 
receding  contact  angle,  any  subsequent  interfacial  motion  would  increase 
the  radius  of  curvature,  resulting  in  less  Laplace  pressure  (which  was 

now  opposing  the  growth  of  the  nucleus)  and  so  the  growth  would  be 

unbounded.  (This  argument  will  be  opposed  in  the  next,  chapter.)  This 
unbounded  growth  would  result  in  a  mostly-vapor-filled  cavity  being 
emitted  from  the  crevice.  This  nucleation  process  must  occur  within  the 
time  frame  cf  an  acoustic  period. 

Using  this  nucleation  criterion  to  define  a  cavitation  threshold, 
the  threshold  can  be  predicted.  Crum  has  modified  previous  results  to 
include  the  effects  of  surface  tension.  He  finds  the  threshold  to  be 

P  -  P  -  P  , 

P  =  (P  -  P  -  P  )  +  — - j - k  { (cosiJO  -l)+(sin<t>)  [  l-(r:  -l)z]  } 

A  L  v  g  o  Y  T 

where : 

P  *  (negative)  acoustic  pressure  PT  »  liquid  pressure 

A  L< 

P  =  vapor  pressure  P  =  gas  pressure 

v  o 

6=  Icos  (a  -S)|  v  =  a  +(3 

A  H 


c  "  constant  “  50  dyn/cm  »  advancing  contact  angle 
6  ■  crevice  half  angle  ct^  »  hysteresis  angle 

This  equation  predicts  the  threshold  as  a  function  of  dissolved  gas 
pressure,  surface  tension,  and  temperature  fairly  well.  More  details  of 
these  predictions  will  be  given  later. 

1-4  Final  Introductory  Remarks 

Each  of  the  models  discussed  above  are  stabilization  mechanisms; 
each  requires  the  presence  of  a  gas  phase  large  enough  to  be  stabilized. 
The  origin  of  the  gas  phase  is  not  of  concern  here.  None  of  the  models 
mentioned  predict  the  subsequent  growth  of  a  nucleus  to  an  observable 
size.  It  is  simply  assumed  that  after  the  bubble  is  nucleated,  it  will 
grow  large  enough  to  be  detected.  This  is,  of  course,  a  huge  assumption 
and  is  not  always  valid.  This  point  will  be  addressed  in  a  subsequent: 
chapter. 

In  this  dissertation,  many  different,  and  sometimes  subtle, 
concepts  are  introduced  and  discussed.  One  way  in  which  these  many 
topics  could  be  handled  would  be  to  discuss  each  in  exhaustive  detail, 
the  result  being  a  much  longer,  and  perhaps  more  complete,  dissertation. 
However,  in  an  effort  to  be  as  concise  as  possible,  this  approach  was 
not  taken.  This  choice  may  lead  to  a  choppy  appearence  of  the  body  of 
the  text.  In  an  effort  to  remedy  this  appearence,  a  brief  description 
of  the  layout  of  this  dissertation  will  be  given. 

Two  types  of  cavitation  will  be  dealt  with.  One  type  is  vaporous 


2b 


(or  transient)  cavitation  in  which  the  cavity  is  mostly  vapor-filled. 
The  other  type  of  cavitation  will  be  termed  "diffusion  cavitation"  and 
it  is  controlled  by  gaseous  diffusion  into  or  out  of  the  bubble.  Three 
different  nucleation  models  will  be  discussed.  Of  the  three,  the 
crevice  model  will  receive  the  majority  of  the  attention.  The  other  two 
models  are  the  ionic  skin  model  and  the  varying-permeability  organic 
skin  model. 

Chapter  2  contains  the  majority  of  the  theoretical  derivations  that 
will  be  presented.  The  concept  of  a  critical  radius  will  be  introduced 
for  both  types  of  cavitation.  This  concept  will  then  be  applied  to  the 
crevice  model  and  expressions  for  both  the  vaporous  and  diffusion 
cavitation  thresholds  will  be  derived.  In  chapter  3,  the  theoretical 
thresholds  predicted  by  the  crevice  model  will  be  compared  to 
experimental  results.  Chapter  '+  is  devoted  to  the  other  models.  The 
plausibility  of  the  ionic  skin  model  will  be  discussed  in  light  of  the 
results  of  measurements  of  the  variation  in  the  acoustic  cavitation 
threshold  with  the  dissolved  ion  concentration.  In  addition,  the 
varying-permeability  model  will  be  applied  to  vaporous  cavitation.  To 
the  author's  knowledge,  this  is  the  first  attempt  at  such  an 
application.  The  final  chapter  will  include  a  summary  of  the  important 
points  presented  in  the  dissertation  as  well  as  tne  conclusions  that  may 
be  drawn  from  them. 
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herivation  of  Cavitatjon  Thresholds 
Using  the  Critical  Radius  Approach  to 
the  Crevice  Model 
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2-1  Introduction 

In  the  previous  chapter,  it  was  noted  that  in  the  current  crevice 
model,  a  vapor  bubble  would  be  nucleated  when  the  liquid-gas  interface 
reached  the  receding  contact  angle.  In  this  dissertation,  it  will  be 
shown  that  this  criterion  is  Indeed  necessary  but  nor.  sufficient.  In 
order  for  a  cavity  to  be  nucleated,  it  is  also  necessary  that  the 
radius  of  curvature  exceed  a  critical  value  determined  liuui  a  stability 
criterion,  the  nature  of  which  depends  upon  the  type  of  cavitation  under 
investigation.  In  view  of  this  requirement,  an  expression  for  the 
liquid  pressure  required  to  produce  such  a  radius  of  curvature  will  be 
derived. 

The  discussion  will  be  divided  into  two  cases:  one  for  ambient 
pressures  less  than  the  vapor  pressure  and  the  other  for  ambient 
pressures  greater  than  the  vapor  pressure.  In  the  first  case  the  cavitv 
formed  will  be  vapor-filled  while  in  the  second  case  the  cavity  will  be 
gas-filled,  growing  by  gaseous  diffusion.  The  results  of  this 
rederivation  will  be  compared  with  those  of  the  previous  crevice  model. 
It  will  be  shown  that  predictions  of  cavitation  thresholds  based  on  this 
revised  crevice  model  agree  closely  with  experimental  results.  In 
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addition,  the  revised  model  wlLl  be  applied  to  situations  previously 
explained  only  in  terms  of  the  varying-permeabilitv  model. 

2-2  Critical  Radius 

The  concept  of  a  critical  radius  is  not  new.  Therefore,  it  Is 

somewhat  surprising  that  this  approach  to  the  crevice  model  has  not  been 

taken  before.  Harvey^  did  discuss  a  critical  radius  for  a  nucleus  to 

grow  bv  gaseous  diffusion.  However,  he  did  not  use  it  in  a 
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quantitative  wav.  Blake  and  Prosperetti  have  discussed  the  critical 
radius  at  which  a  vapor  bubble  is  stable  against  growth  or  collapse. 
Building  on  these  concepts,  expressions  will  be  derived  for  the  applied 
pressure,  called  the  cavitation  threshold,  reauired  for  a  bubble  to  be 
nucleated  from  a  crevice.  (Wore:  cue  applied  piessuie  mav  assume 
negative  values.)  Since  the  spherical  geometry  of  a  free  bubble  lends 
itself  to  easier  manipulation  than  the  conical  geometry  of  a  crevice, 
the  critical  radii  will  be  derived  for  a  free  bubble.  In  general,  the 
stability  conditions  hold  lor  either  case;  exceptions  will  be  noted  when 
necessary. 

2-2.1  Critical  Radius  for  Growth  by  Caseous  Diffusion 

Consider  a  gas  filled  bubble  of  radius  R  surrounded  by  a  liquid 
which  exerts  a  pressure  on  the  bubble.  The  liauid-gas  surface 
tension  is  .  The  gas  pressure  inside  the  bubble  is  P  ;  the  vapor 

g 

pressure,  p  ,  is  assumed  negligible  (see  figure  6).  Assume  the  bubble 
is  in  equilibrium  so  that  the  pressure  inside  (P  )  equals  the  liquid 

g 

pressure  plus  the  Laplace  pressure  due  to  the  surface  tension: 


Figure  6.  Illustration  of  a  free  bubble  surrounded 
by  a  liquid.  The  external  pressure  is  the  sum  of 
the  liquid  and  Laplace  pressures  and  the  internal 
pressure  is  the  sum  of  the  gas  and  vapor  pressure. 
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F  =  PT  +  .  ( 1 ) 

g  L  R 

Simultaneously,  the  svstera  roust  be  in  diffusion  eaui 1 ibr iuro.  This 

requires 

c.  =  a  c.  (2) 

L  l 

where  c  and  c  refer  to  the  concentrations  of  gas  dissolved  in  the 
liquid  and  present  inside  the  bubble,  respectively;  "a"  is  the 
solubiHty  coefficient  of  the  gas  in  the  liquid.  Tf  the  gas  is  ideal 
then  (2)  becomes 


where  B  is  the  universal  gas  constant  and  T  is  the  absolute  temperature. 

Equation  (3)  is  a  form  of  Henrv's  Law.  At  equilibrium  P  mav  be 

E 

replaced  with  t,  called  the  gas  tension,  which  is  a  measure  of  the 

partial  pressure  of  the  gas  dissolved  in  the  liquid.  At  equilibrium 

P„  =  Pt  +  41  =  1  <*) 

g  L  R 


<“  -  V  -  'TT 

The  critical  (equilibrium)  radius  is,  therefore, 


The  subscript  "eg"  refers  to  the  critical  radius  for  a  bubble  to  grow  by 
gaseous  diffusion.  Notice  that  the  radius  of  curvature  is  positive 
only  in  supersaturated  solutions.  Since  a  negative  radius  is 
meaningless  for  free  bubbles,  stabilitv  can  occur  only  in  supersaturated 
solutions. 


i . . 


31 


2-2.2  Critical  Radius  for  Crowth  in  the  Absence  of  Diffusion 

The  discussion  presented  in  this  section  is  similar  to  that  of  the 
previous  one,  except  now  the  effect  of  vapor  pressure  is  included. 
Again,  at  equilibrium  the  following  condition  holds 


Pv+Pg“PL+-X*  (7) 
The  condition  given  by  equation  (3)  is  dropped  and  the  amount  of  gas 
inside  the  bubble  is  fixed.  (In  other  words,  no  diffusion  of  gas  takes 
place.  However,  the  amount  of  vapor  is  not  fixed.)  If  the  gas  is  ideal 


then  rearranging  (7)  gives 


p  +2l,  3S»I+  p 

L  R  4irRJ  v 


where  n  ic  the  number  cf  rue i s s  of  the  gas  inside  the  bubble e  To  examine 
the  stability  consider  the  effect  of  a  small  increase  in  R.  For  the 
bubble  not  to  grow  spontaneously,  the  effective  external  pressure,  which 
is  given  by  the  left-hand  side  of  equation  (8),  must  increase  more  than 
the  internal  pressure,  given  by  the  right-hand  side  of  equation  (8). 
This  will  always  be  true  if  _>  P  ,  because  the  gas  pressure  decreases 
more  with  an  increase  in  R  than  does  the  Laplace  pressure.  If,  however. 


P,  <  P  then  the  situation  is  different. 

L  v 

The  stability  criterion  can  be  expressed,  in  general,  as 

—  (p  +  iLl)  >  JL  (  +  p  )  .  (9) 

9R  V  I.  R  '  9R  '  4uR3  V 

Performing  the  differentiation  and  solving  for  R  gives 

R  <  irp^hnr  •  (10) 

'  V  L 

Tn  other  words,  when  >  P^,  the  bubble  will  be  stable  against 


t 


spontaneous  growth  only  if  R  <  R  where 

Rcv  -  7(Pv4-  PL)  •  “*> 

Taking  the  other  view,  a  bubble  will  grow  in  a  spontaneous  unbounded 

manner  when  the  liquid  pressure  is  less  than  the  vapor  pressure  only  if 

R  >_  Rcv*  A  bubble  that  grows  in  this  way  will  be  mostly  vapor-f i lied  * 

Therefore,  the  subscript  "cv"  refers  to  the  critical  radius  for  a  vapor 

bubble  to  grow  from  a  gas-filled  nucleus.  The  reader  is  referred  to 
38 

Prosperetti  for  more  details. 


2-3  Threshold  Derivations 


In  the  preceding  section,  the  conditions  for  the  nucleation  of  both 
gaseous  and  vaporous  hubbies  were  developed.  Using  them,  expressions 


V.v'.'vVi 


will  now  be  derived  for  the  applied  pressures  required  for  nucleation 
from  gas-filled  crevices.  These  are  referred  to  as  cavitation 
thresholds.  It  should  be  noted  that  the  applied  pressure,  P^,  is  equal 
to  the  liquid  pressure  minus  the  ambient  static  pressure;  i.e., 
P  *  P  -  P  .  (P  experl  illy  is  the  applied  acoustic  pressure.)  The 

n  L  O  A 

discussion  which  follows  will  be  divided  into  two  different  cases:  one 
for  PT  <  P  and  the  other  for  P.  >  P  .  This  division  is  a  natural  one; 
the  difference  in  the  restrictions  on  is  in  a  sense  the  difference 
between  acoustic  cavitation  and  diffusion  cavitation.  Acoustic 
cavitation  is  produced  by  negative  applied  pressures  and  its  nucleation 
usually  involves  the  growth  of  vapor-filled  bubbles,  although  not 
always.  On  the  other  hand,  diffusion  cavitation  is  usually  produced  by 


ambient  pressures  greater  than  the  vapor  pressure  and  gas  diffusion  is 
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the  predominant  mechanism  in  its  nueleation.  Diffusion  cavitation 
should  not  be  confused  with  gaseous  (stable)  acoustic  cavitation;  the 
two  occur  in  different  regions  of  ambient  pressures.  Acoustic 
cavitation  will  be  considered  first. 


2-3.1  Acoustic  Cavitation  Nueleation 

Assume  that  a  deep  crevice,  partially  filled  with  gas,  exists  in  a 
hydrophobic  solid  surrounded  by  a  gas-saturated  liquid.  The  hydrostatic 
pressure  due  to  the  liquid  is  negligible  and  it  is  exposed  to  an 
atmosphere  at  pressure  Pq  (refer  to  figure  7).  The  liquid-gas 
interface  contacts  the  crevice  wall  at  a  height  z  above  the  apex.  The 
crevice  angle  is  22  s  The  system  is  in  equilibrium,  so 


P  «  t  -  P  . 
o  g 


(12) 


For  temperatures  far  below  the  boiling  point,  P^  «  Pq  and  thus  the 
interface  is  essentially  flat. 

The  liquid  is  now  degassed  so  that  t  <  Pq.  Gas  will  diffuse  out  of 

the  nucleus  causing  the  interface  to  bow  inward,  toward  the  apex,  as 

shown  by  the  dashed  line  (b)  in  figure  7.  At  equilibrium  P  ■  t  and 

g 


p  +  P  +  11 
g  V  R 


(13) 


The  new  equilibrium  position  of  the  nucleus  is  shown  by  the  solid  line 
(c)  in  figure  7.  At  this  position  the  angle  between  the  interface  and 


the  wall  is  the  advancing  contact  angle. 

At  this  point,  an  acoustic  field  of  pressure  amplitude  is  applied 
such  that  the  liquid  pressure  exerted  on  the  interface  is  PQ  +  P^. 


( 
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Assume  is  such  that  -  P^.  During  the  negative  half-cycle, 

the  interface  bows  outward  and  a  new  equilibrium  position  is 

established.  This  process  will  be  isothermal  if  the  thermal  diffusion 
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length,  L,  is  large  compared  to  crevice  dimensions.  L  is  given  by“ 

L  -  <— )** 
pc  u/ 

P 

where  k  is  the  thermal  conductivity  of  the  gas,  p  is  its  density,  c^  is 
the  specific  heat  at  constant  pressure  per  unit  mass,  and  w  =»  2nf  where 
f  is  the  frequency  of  the  acoustic  field.  For  most  gases  at  normal 
pressures  and  temperatures  (P  s  1  bar  and  T  5  200-500K),  k  is  on  the 
order  of  5x10  ^  kcal  m  *s  *K  * ;  the  density  of  air  is  1.1  kg  m  3 ;  c 

P 

is  0.25  kcal  kg  Therefore, 


L  s  2,4  v  10-3  f~i/2  m  .  f  1 5 ) 

Notice  that  L  is  inversely  proportional  to  the  square  root  of  the 
acoustic  frequency,  but  even  at  megahertz  frequencies,  L  »  2.4  y m.  This 
is  about  an  order  of  magnitude  larger  than  the  crevice  sizes  assumed 
here.  Therefore,  the  isothermal  assumption  is  a  good  one. 

A  bubble  will  be  nucleated  flora  this  new  position  only  if  the  radius 
of  curvature  of  the  interface  is  larger  than  the  critical  radius  R 

cv 

The  value  of  P^  necessary  for  this  condition  to  be  met  is  found  as 
follows. 

For  isothermal  processes, 


P  V  =  P  V 
11  2  2 

where  P^  and  are  the  internal  gas  pressure  and  nuclear 
initial  degassed  state  whereas  and  are  the  interior 


and  nuclear  volume  at  P  ,*  P  -  P„.  V, 

L  o  A  1 


(16) 

volume  in  the 
gas  pressure 


is  given  by 


1 


cone 


cap 


(see  figure  8).  After  some  geometric  and  algebraic  manipulation 

tt  =  -ir-a3  (cut3-n  )  (18) 

where 

n 

i 


and 

6^  =  | cos(ax  -  8) |  • 

(Throughout  this  dissertation  the  subscripts  associated  with 
refer  to  the  subscript  associated  with  u.)  The  final  volume 
similarly  given  by 

V2  “  I™* (cotg  +  n2)  * 

Tlie  interior  pressures  are  given  by 


(20) 

n  and  6 

V  is 


(21) 


P  -  P  (22) 

1  g 


and 


2  Y 


P  =  P  -  P  -  P  +  -5-  . 

2  o  A  v  R 


(23) 


Rearranging  (16)  gives 


(24) 


or 


P  -  P  -  P  +  ~  (cot8-n  ) 

o  A  v  R?  _  3  i _ i 


-=-na  (cotf3+n  ) 
3  2  2 


(25) 


Using  the  geometrical  relationship  that  a  =  OR,  (25)  becomes  (after  some 
algebra) 


R 


3 


i~  cot6-n  -i 
_ 1. 

-  cot  8+0  -1 


=  0  . 


(26) 
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At  the  initial  state,  P.  =  P  »  t.  From  equation  (13)  it  is  seen  that 

1  g 


2y 


or 


Ri  ~  p  -  P  -  P 
o  g  v 

2y6 


a  = 


i  P  -  P  -  P  ' 
o  g  v 


Using  this  expression  and  letting 

R  «  R 

gives 


±L 


2  cv  3<Pv+Pa-Po) 


(27) 


(28) 


(29) 


P  =  P  -  P  +  P  .  , 
A  o  v  ,  V2  g  L! 


/  P  -  P 
I  o 


H 


r  cot3+n 


[  cot^-n 


(30) 


In  general,  In  the  degassed  state,  the  contact  angle  is  the  advancing 


contact  angle, 


In  the  final  state  the  contact  angle  equals  ths 


receding  contact  angle,  a^-a^.  Equation  (30)  new  becomes 

[cotB+n. 

PA  Po  Pv  +  i/2  Pg  l_[  P 


-  1 


J  'cos(aR  -B) 


cos (a  -8) 
A 


[cotS-n  jj 


rr 

.  (31) 


This  relationship  gives  the  threshold  for  vaporous  cavitation.  Although 

a  more  detailed  comparison  will  be  made  in  the  next  chapter,  a  few  brief 

comparisons  to  the  old  crevice  model  should  be  made  here.  Apfel's 
34 

threshold  is 


P  =  P  -  P  -  P  +  (P  -  P  -  P  ) 
A  o  v  g  o  v  g' 


This  equation  gives  a  linear  dependence  on  gas 

not.  The  I  cos  (a  -S)/cos  (a  ~3)|  term  appears 
R  A 

is  different. 


cos (a  -B) 
K _ 

cos (a  -B) 

A 


(32) 


content  whereas  (31)  does 
in  both  but  the  dependence 


i. 
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It  was  assumed  throughout  the  derivation  that  the  interface  reaches 
the  receding  contact  angle.  For  small  values  of  P  this  may  not  be 

ri 

true.  In  order  to  find  the  value  of  necessary  for  the  interface  to 

reach  a  ,  again  assume  isothermal  conditions..  The  initial  state  is  the 
K 

degassed  stated  mentioned  above  with 

P.  ”  (33) 

1  g 


v,  "  (cotg  -n . )  • 

1  3  1  A 

The  final  state,  2,  is  such  that  the  interface  just  reaches  o  so 

R 

2y 

p  *  p  „  p  -  p  +  --L 

2  o  v  A  R2 
V  «  vira3  (cotB+r>_)  . 


Recall  that 


a  *  5  R 


P  “P-P-P,  +  - 

2  o  v  A  a 


Using  (16), (28),  and  (37)  it  can  be  shown  that 

6  cot &-n 

V  V  PA  +  £-  (V  V  V  *  P8 


cotB-n, 


cos (a  -B) 

R  (40) 


P  -  p  -  p  - -  +  (F  -  P  -  P  ^  - - - 

o  v  g  (cot0+nD  J  o  v  gJ  cos (a  -&) 

K  '  A 


The  asterisk  refers  to  the  fact  that  this  is  the  critical  value  of  P. 

e 

below  which  nucleation  will  not  occur.  (In  this  dissertation  ai 


asterisk  will  always  denote  some  type  of  "critical"  pressure.)  Why  is 


40 


this  a  critical  value?  Assume  the  outwardly  bowed  interface  has  a 

radius  of  curvature  greater  than  R  ,  vet  it  has  not  reached  the 

cv 

receding  contact  angle.  The  nucleus  is  unstable  and  so  the  radius  of 

curvature  will  not  spontaneously  decrease.  However,  this  would  be 

required  In  order  for  the  contact  angle  to  decrease  to  Since  the 

interface  will  not  move  before  it  reaches  the  receding  contact  angle, 

there  is  no  nucleation.  It  should  be  noted  that  equation  (40)  is  almost 

identical  to  equation  (32),  Apfel's  threshold  equation.  The  only 

difference  is  that  in  deriving  equation  (32)  it  was  assumed  at  gas 

diffusion  can  keep  up  with  the  motion  of  the  interface.  This  assumption 

was  not  made  in  deriving  equation  (40),  hence  the  additional  factor 

associated  with  P  .  If  gas  diffusion  can  keep  up  with  the  motion  of  the 
8 

interface  then  that  factor  equals  one. 

The  implications  of  these  thresholds  will  be  discussed  in  detail  in 
the  next  chapter.  We  turn  now  to  the  threshold  for  the  nucleation  of  a 
bubble  by  gaseous  diffusion. 

2-3.2  Diffusion  Cavitation  Nucleation 

Once  again  assume  a  deep  crevice,  partially  filled  with  gas,  exists 

in  a  hydrophobic  solid  surrounded  by  a  gas-saturated  liquid  that  exerts 

a  static,  ambient,  pressure  on  the  interface.  Tn  equilibrium,  the 

pressure  balance  equation  (12)  holds;  the  interface  is  essentially  flat 

and  a  height  zq  above  the  apex.  Assume  the  ambient  pressure  is 

increased  from  P  to  P  .  This  change  in  pressure  appears  at  the 
o  m 

interface  as  an  increase  in  the  liquid  pressure  causing  it  to  bow  inward 


i 


(see  figure  9).  The  interface  will  advance  toward  the  apex  if  it 

reaches  the  advancing  contact  angle.  The  pressure  at  which  this  occurs 

is  found  by  assuming  an  isothermal  compression  from  V.  -  V  to 

°  j.  cone 

V  -  V  -  V  .  That  is;  from 
f  cone  cap 

V  =  -^a3cote 


V,  =  -^ira3  (cot0-n  )  . 
f  3  A 


The  internal  pressure  also  changes  from  P^  =•  ■  P^  (neglecting  P^)  to 


2  y  2vSa 

?c  =  p  -  ±-L  =  p _ A 

f  m  R  m  a 


Using  P.V  "  P -V  ,  we  have 
l  i  f  f 


1~j(1  ~  uA  tanfs) 


Let  P  .  =  P  -  P  (see  reference  26).  Substitution  into  equation  (44) 
crush  m  o 

gives 


p*  ,  P  -i 

crush  o  [_  1-n  tanS  J 


where  the  asterisk  indicates  that  this  is  the  critical  value  of  P 


crush 


Let  the  crevice  mouth  have  a  radius  A  and  define  the  ratio  a/A  as  h. 
This  non-dimensional  parameter  is  a  measure  of  the  degree  of  fullness  of 
the  crevice.  From  the  conical  geometry  it  is  seen  that 


a  £  _  ^ 


where  z^  is  the  depth  of  the  crevice.  Using  this  relationship  the 
expression  for  the  critical  crushing  pressure  becomes 


the  compression  phase  c-r  a  o 
process  (diffusion  cavitatio 


A3 


crush 


nAtang 

1-0 . tang 
A 


2y6a 

hA 


( A7  ) 


For  crushing  pressures  less  than  this  value,  the  interface  does  not 

reach  the  advancing  contact  angle  and  therefore  cannot  move. 

* 

If  the  crushing  pressures  exceed  P  c  the  interface  advances 

until 


P 

m 


(48) 


The  prime  indicates  that  the  gas  pressure  is  no  longer  equal  to  c,  but 

is  greater.  Sim  a  the  liquid  is  assumed  incompressible ,  the  change  in 

liquid  pressure  does  not  cause  an  immediate  change  in  t ,  as  it  dees  in 

P  .  The  change  in  gas  tension  occurs  at  a  rate  controlled  by  the 
S 

diffusion  of  gas  from  the  atmosphere  into  the  liquid.  Thus,  for  a  time, 

the  nucleus  is  supersaturated.  Gas  diffuses  out  of  the  nucleus  until, 

once  again,  P  =  t.  At  this  point,  the  interfacial  radius  of  curvature 
8 

is  less  than  R.  because  P  <  P  Therefore,  slice  the  interface  is 

l  g  g 

already  at  ,  this  reduction  in  radius  implies  that  the  diffusion 
causes  the  nucleus  to  contract;  that  is,  h  decreases.  At  this 
equilibrium  position 


P  =  P  + 
m  g 


2 


(49) 


where  P  <  P  ’  arid  R.  <  R,  . 
g  g  2  1 

If,  at  this  point,  the  ambient  pressure  were  reduced  back  to  P  , 

o 

the  nucleus  would  expand,  but  not  to  its  original  volume  since  there  are 
fewer  gas  molecules  within  the  nucleus.  This  net  reduction  in  the 
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90 
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volume  of  the  nucleus  would  result  in  a  higher  cavitation  threshold  than 

that  of  an  identical,  non-pressurized  nucleus*  Such  an  increase  in 
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threshold  is  well  documented.  *  If  the  liquid  pressure  is 

maintained  at  long  enough,  however,  the  dissolved  gas  pressure  will 

eventually  reach  P  ;  that  is,  t  =  P  .  This  results  in  the  diffusion  of 
ni  mm 

gas  back  into  the  nucleus  until  P  =  P  ,  at  which  point  the  interface  is 

g  in  * 

essentially  flat  and  a  height  above  the  apex. 

Suppose  now  the  atmospheric,  and  therefore  the  liquid,  pressure 
were  reduced  to  a  final  pressure  P^.  The  internal  gas  pressure  of  the 
nucleus  will  be  greater  than  the  external  ambient  pressure;  the 
interface  bows  outward.  In  order  for  the  interface  to  reach  the 
receding  contact  angle  the  volume  must  change  from 


V  =  -Wa3  cotB 
1  3  2 


(50) 


to 


V  =  -^-ira3  (cotg+nj 

I  d  2  K 


(51) 


while  the  gas  pressure  inside  the  nucleus  chanses  from  P  to 

2yd 


p  +  «  p  +  R 

'  “  f  h  A  ’ 

2 


f  '  R 

Assuming  isothermal  conditions,  we  have 

2yd„ 


p  = 
m 


P  +  - — - 

f  h  A 


(1+u  tanb)  . 


(52) 


(53) 


The  extent  to  which  the  nucleus  is  supersaturated  can  be  characterized 
by  defining  the  supersaturation  pressure  ?gc  =  -  P^.  Using  this 


expression  and  noting  that 


F=P  +  P  -  P 

f  o  crush  ss 


(54) 


< 
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equation  (53)  becomes 


P  =  (P  +  P  .  ) 
ss  o  crush 


r  r,Rtan6 

l+nr  tan£ 
K 


H' 


h  A 

?. 


(55) 


Notice  the  similarity  between  this  equation  and  equation  (47).  The 
difference  in  sign  in  the  denominator  arises  from  the  fact  chat  in  one 
instance  the  calotte  adds  to  the  volume  of  the  nucleus,  while  in  the 
other  it  subtracts  from  it. 

For  supersaturation  pressures  greater  than  P  ,  the  interface 

reaches  the  receding  contact  angle  ar.d  moves  outward.  If  the  change  in 

volume  is  such  that  the  radius  of  curvature  of  the  interface  is  larger 

than  R  ,  equation  (6),  the  nucleus  will  grow  by  gas  diffusion.  la 
eg 

general ,  the  final  volume  (51)  is  gJ ven  hv 


Vf  “  7*a;<cotf5+nR) 


where  >  a2* 


Recall  that 


In  this  case  (53)  becomes 


P  - 
m 


P  +  h 

f  R 


ca ) 


a 

a 

2J 

0+nR  tang) 


eg 


2y 

t-P. 


(56) 


(57) 


(53) 


In  this  case,  t  *>  P  end  P.  *  P-.  Therefore, 
m  l  r 

R  = 
eg  P 
°  ss 

By  definition.  P  -  P,  +  P  .  Therefore,  equation  (57)  reduces  to 
J  ’  m  f  ss 

>  3 


fV 

(32j 


(1+n  tanS)  »  1 


(59) 


(60) 


Also  note  that  a,  -  6  R„  where  R-  is  defined  in  (49).  Equation  (49)  can 
z  A  2  2 


be  rewritten  as 
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- -  K  P  -  P  =  P 

R-  mg  crush*  (61) 

2 

Substituting  (61)  into  (60)  and  recalling  that  “ipiRCg»  one  finds 


ss  =  Pcrush  lT1|(1+Van6)3 


lcos(a_-g)  A. 

P  ■  P  | - — -  (1  +  n  tanSf 3  (63) 

ss  crush  |cos(a  -B)  ',1  +  nRt  * 

A 

Tills  is  the  threshold  for  the  growth  of  a  nucleus  by  gas  diffusion 

*  * 

provided  that  P  >  P  and  P  .  >  P  ,  . 

ss  ss  crush  crush 

Before  getting  into  the  analysis  of  the  model  equations  presented 
here,  it  is  perhaps  worthwhile  to  dwell  a  bit  longer  on  the  concept  of  a 


critical  crushing  oressure.  P  ,  is  given  by  equation  (47): 
~  "  crusn 


P*  -  P  r  vanB  1 1 

crush  ol  1-n  tang 

'-A  --1 


A 

hA  • 


Solving  this  for  the  radius  of  the  crevice  mouth  gives 


2yo  _  *  n.tanjj  |  i  * 

=  - _  p  -  p - - 

h  I  crush  o  1-r.  tang  I 
u  A  1  J 


p  tang 


n  a 

Writing  P  ,  as  P  -  P  »  (65)  becomes 
°  crush  m  o 

2^a  r*  f 

A  - - -  P  —  P  i 


■*  n  1  1 

F  -  P  - - - 1 

m  o  1-n^tang 


From  equation  (44)  one  sees  that  t.he  second  term  inside  the  brackets  is 

simply  equal  to  the  final  gas  pressure,  F' ,  inBide  the  nucleus.  Using 

g 

this  ( 66)  becomes 

7  |cos(a  -3)1 

A  =  Ll  J___A__L  (67) 

h  (P*-  P'  )  • 

m  g 


Apfel  '  defines  a  critical  crevice  radius  for  acoustic  cavitation  as 


l 
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*  2-y  |  cos(a^-6) 

A.  =  /  p  _p  ^  ( 68 ) 

L  g 

* 

(neglecting  vapor  pressure).  Apfel's  P.  and  P  are  analogous  to  P  and 

L  g  m 

* 

in  equation  (67).  Therefore,  the  concept  of  P  crugh  in  diffusion 

* 

cavitation  is  analogous  to  Apfel's  concept  of  A  in  acoustic  cavitation, 
the  only  difference  being  that  Apfel  considers  only  full  crevices 
(h  -  1). 


c 


c 


Chapter  3 

Analysis  of  the  Critical  Radius 
Approach  to  the  Crevice  Model 


3*1  Introduction 

In  this  chapter  the  critical  radius  approach  to  the  crevice  model 
will  be  analyzed  on  the  basis  of  its  ability  to  predict  cavitation 
thresholds.  The  discussion  will  include  both  acoustic  and  diffusion 
cavitation  and  only  the  crevice  model  will  be  considered.  In  a  later 
chapter  other  stabilization  and  nucleation  models  will  be  analyzed. 


3-2  Application  of  the  Crevice  Model  to  Acoustic  Cavitation 

The  acoustic  cavitation  threshold  of  water  will  be  predicted  as  a 
function  of  dissolved  gas  pressure,  surface  tension,  and  temperature. 
Unless  otherwise  stated,  all  theoretical  predictions  will  be  based  on 
the  critical  radius  criterion. 


3-2.1  The  Effect  of  Dissolved  Gas  Pressure 

Recall  equations  (31)  and  (AO)  from  the  previous  chanter: 


,  P  -  P 

2  o  v 

?  -  P  -  P  +  ~  P  — r - 1 

A  o  v  3H  g  LI  P 

and  g 


J  cos(o  -B)  3  cot0+n 

iv.  K 


[  cos(a^  — B) I  (cotB-nA 


*  cotB-n 

P  -  p  -  p  -  p  - — 

A  o  v  g  ^cotB+nR 


c.os(ot  -B) 

+  (P  -  P  -  P  )  - >— - -T- 

o  V  g  cos (a  -B) 


These  equations  are  the  expressions  for  the  acoustic  pressures  necessary 

for  the  interface  to  reach  the  critical  radius  (R  )  and  the  receding 
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contact  angle,  respectively.  Equation  (AO)  is  essentially  the  same  as 


A8 


the  threshold  predicted  by  the  previous  formulation  of  the  crevice 

model. 

The  dependence  of  these  thresholds  on  the  dissolved  gas  pressure  is 
explicitly  stated.  (We  shall  see  later  that  the  dependences  on  surface 
tension  and  temperature  are  not  so  explicit.)  Recall  that  P^  is 

equivalent  to  the  partial  pressure  that  the  gas  inside  the  nucleus 

exerts  on  the  interface  in  both  the  saturated  and  ur.dersaturated 
(degassed)  states.  In  the  absence  of  gas  diffusion,  the  partial 

pressure  of  the  gas  in  an  expanding  nucleus  is  less  than  P^. 

Figures  10a  and  10b  show  the  cavitation  threshold  as  a  function  of 
dissolved  gas  content  for  two  representative  values  of  surface  tension. 
These  graphs  illustrate  several  important  points.  First,  the  two 
nucleation  criteria  are  equivalent  over  a  short  range  of  gas  pressures. 
This  equivalence  is  difficult  to  show  analytically;  however, 
numerically  it  is  easily  demonstrated.  Hence  for  threshold  measurements 
made  in  these  regions,  the  different  criteria  give  equivalent 
predictions.  The  theories  diverge  on  either  side  of  this  equivalent 
region.  At  no  value  of  gas  pressure  is  the  threshold  predicted  by  the 
critical  radius  criterion  (P  )  less  than  that  predicted  by  the  receding 

■k  it 

contact  angle  criterion  (P  ).  In  other  words,  while  P  is  a  necessary 

A  A 

condition  for  nucleation,  P^  is  both  a  necessary  and  sufficient 
condition  for  nucleation. 

At  low  gas  pressures,  the  receding  contact  angle  criterion  gives  a 
limiting  threshold  value  of  between  about  6  and  12  bars  for  the  ranges 


of  surface  tensions  discussed  here.  The  exact  values  are  not  important. 


What  is  important  is  that  the  values  are  relatively  small  compared  to 
the  threshold  for  homogeneous  nucleation  as  discussed  in  the  first 
chapter.  This  relatively  low  prediction  is  not  a  realistic  one.  It  has 
been  assumed  that  the  nuclei  are  gas-filled;  the  fact  that  pre 
pressurization  and  variation  of  gas  pressure  have  experimentally  been 
shown  to  have  an  effect  on  thresholds^  verifies  this  assumption  of  gas- 
filled  nuclei.  It  follows,  therefore,  that  if  all  of  the  gas  is  removed 
from  the  water,  that  is  if  the  gas  pressure  is  decreased  to  zero,  the 
nuclei  would  in  effect  be  removed.  Hence,  the  threshold  for  completely 
degassed  systems  should  approach  homogeneous  nucleation  thresholds. 
This  expected  behavior  is  predicted  by  the  critical  radius  criterion. 

Figures  11a  and  lib  are  two  graphs  of  the  cavitation  threshold 

versus  gas  pressure.  Plotted  are  the  theoretically  predicted  thresholds 
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as  well  as  thresholds  measured  by  Strasberg  ,  Connolly  and  Fox  ,  and 
1 8 

Galloway  .  No  attempt  has  been  made  to  normalize  any  of  the  data, 
although  the  calibrations  of  the  various  experimental  systems  were  based 
on  different  standards.  Three  comparisons  are  evident.  If  the  "lone" 
data  point  at  a  gas  pressure  of  about  0.1  bar  were  eliminated,  there 
would  appear  to  be  a  linear  dependence  on  gas  pressure.  However,  the 
"lone"  data  point  is  not  actually  "lone".  Figure  lib  is  an  expanded 
version  of  figure  11a  which  includes  Galloway's  data  for  gas  pressures 
below  0.1  bar.  (These  are  the  only  data  available  for  these  extreme 
values  of  gas  pressure.)  Emphasis  should  not  be  placed  on  the  absolute 
values  of  his  data;  thresholds  of  those  magnitudes  are  uncommon.  The 
trend  of  the  data,  however,  is  important  and  in  good  agreement  with  the 


re  11a.  Graph  of  threshold  versus  gas  pressure 
Strasberg's  (U) ,  Galloway's  (0) »  and  Connolly  a 


(JDq)  anOHS3dH 


ure  lib.  Graph  of  threshold  versus  gas  pressure  for 
pressures.  The  data  are  Galloway's  and  the  solid  11 
the  prediction  of  the  critical  radius  approach  to  the 
vice  model. 


critical  radius  prediction 


It  should  be  pointed  out  that  the  limiting 


value  of  the  threshold  as  predicted  by  the  receding  contact  angle 

criterion  differs  from  the  experimental  and  theoretical  values  given  in 

figure  ] lb  by  well  over  an  order  of  magnitude. 

As  shown  in  figure  10,  the  receding  contact  angle  criterion 

predicts  a  threshold  cf  zero  bars  for  near-saturated  solutions.  On  the 

other  hand,  the  critical  radius  criterion  predicts  a  threshold  near  1 

bar  (actually  P  — P  ,  which  is  approximately  1  bar  for  low  temperatures), 
o  v 

The  difference  is  inherent  in  the  definition  of  R  (given  by  equation 
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11  of  the  previous  chapter)  in  that  the  expression  for  R  is  valid 
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onlv  for  total  external  pressures  (P  -P.)  less  than  P  •  There  is  no 

o  A  v 

k 

such  restriction  on  P  .  This  difference,  although  forced,  should  be 
experimentally  verifiable. 

Referring  to  figure  Ha,  it  is  seen  that  the  experimental 
thresholds  tend  to  be  about  2  bars  for  saturated  solutions.  Although 

A 

this  is  in  closer  agreement  with  P.  than  P  the  discrepancy  is 

A  A 

troublesome.  This  discrepancy  may  result  from  the  fact  that  the 
theoretical  model  does  not  strictly  reproduce  the  experimental 
procedure.  That  is,  the  lines  on  the  graphs  represent  the  acoustic 
pressure  amplitude  necessary  for  the  nucleus  to  become  mechanically 
unstable.  The  experimental  thresholds  are  measured  values  of  the 
acoustic  pressure  amplitude  within  the  liquid  at  the  time  a  bubble  was 
detected.  It  should  not  be  assumed  that  these  two  pressure  amplitudes 
are  always  the  same.  The  process  by  which  a  nucleated  bubble  (or 
nucleate)  becomes  a  detected  bubble  is  indeed  complex,  ns  alluded  to  at 
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the  close  of  the  first  chapter.  One  factor  which  can  complicate  the 
evolution  of  a  nucleate  is  gas  diffusion  into  the  expanding  bubble. 
This  factor  is  ignored  in  the  theoretical  threshold  predictions. 
However,  the  effects  of  gas.  dif fusion  should  diminish  as  the  gas  content 
of  the  system  decreases.  This  conclusion  is  consistent  with  the  results 
shown  in  figure  11a;  there  is  better  agreement  between  theory  and 
experiment  at  gas  pressures  below  about  0.5  bar. 

A  somewhat  quantitative  explanation  of  the  discrepancy  can  be  made 

35 

with  the  use  of  Apfel's  work  on  acoustic  cavitation  prediction.  Apfel 

defined  thresholds  for  various  types  of  cavitation  and  uses  these 

definitions  to  predict  which  of  the  types  of  cavitation  will  occur  under 

a  given  set  of  experimental  conditions.  The  thresholds  reproduced  in 

figure  11  were  described  by  the  experimenters  as  corresponding  to 

"transient"  cavitation  events.  Tf  an  order  of  magnitude  estimate  of  the 

initial  radius  of  a  nucleate  is  made,  Apfel's  criteria  for  transient 

cavitation  can  be  applied  to  find  the  pressure  amplitude  necessary  for  a 

nucleate  to  evolve  into  a  detected  bubble.  The  term  "detected  bubble" 

refers  to  a  bubble  from  which  a  vaporous  cavity  can  form  in  the 

presence  of  a  sound  field  of  sufficient  amplitude. 

Up  to  this  point  little  has  been  said  concerning  the  physical  size 
92 

of  the  crevice.  Apfel  estimated  the  radius  of  a  crevice  mouth  to  be 

about  0.5  um.  Assume  that  the  crevice  has  a  depth  equal  to  its  mouth 

-19  3 

diameter.  This  results  in  a  crevice  volume  of  about  3x10  m  .  Let  us 


further  assume  that  the  nucleation  process  deactivates  the  crevice.  By 
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this  it  is  meant  that  when  the  critical  radius  is  reached,  the  expanding 

vapor  cavity  will  eventually  detach  itself  from  the  solid,  removing  all 

of  the  gas  from  the  crevice  in  the  process.  Thus ,  shortly  after 

nucleation,  a  mostly-vapor-filled  cavity  exists  which  contains  an  amount 

of  gas  such  that  under  average  conditions  (zero  acoustic  pressure)  would 

-19  3 

have  a  volume  of  3x10  m  *  Therefore  in  order  to  predict  the  future 
of  this  nucleate,  we  shall  examine  the  evolution  of  a  spherical  bubble 
having  that  volume,  or  equivalently,  a  bubble  having  an  initial 
equilibrium  radius  of  about  0.4  pm. 

Apfel's  analysis  indicates  that  a  bubble  small  compared  to  its 
resonance  radius  will  become  transient  when  the  acoustic  pressure 
amplitude  exceeds  the  Blake  threshold.  The  Blake  threshold  is  the  free 
bubble  analogy  to  the  critical  radius  approach  to  the  crevice  model.  It 


is  derived  by  combining  the  value  of  R  with  the  pressure  balance 
equation  for  a  f'-ee  bubble.  The  Biake  threshold  is 

4  3XB  l4~j  (69) 

PB  =  Po  1  +  9  XB  (4 ( 1+X  ) j  J 
whe  re  ° 


9 


P  R  (70) 

o  o 

Here,  Rq  is  the  equilibrium  size  of  the  bubble.  Using  P  ~1  bar,  y=70 
dyn/cm,  and  R  =0.4  pm,  the  calculated  value  of  P_  is  2.2  bar.  Figure 

O  D 

lla  has  been  reproduced  in  figure  12  with  the  exception  that  Pn  is  shown 

a 

* 

instead  of  P  .  The  lines  intersect  at  a  gas  pressure  of  about  0.6  bar. 

c\ 


Below  this  value  PA  is  larger  and  the  fit  with  the  data  is  good  in  this 
region.  However,  above  0.6  bar  the  prediction  of  diverges  from 


.  '.1 


experiment  and  it  is  seen  that  P„  is  a  much  better  predictor.  It  should 
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be  pointed  out  that  the  straight  solid  line  assumes  a  constant  value  of 

R  . 

o 

This  analysis  of  thresholds  near  saturation  does  not  prove  that  the 
cavitation  nucleus  is  in  fact  a  gas-filled  crevice.  It  is,  however,  not 
inconsistent  with  the  hypothesis.  It  is  seen  that  a  gas-filled  crevice 
could  give  ri  e  to  a  nucleate  which  would  later  grow  as  a  free  bubble  to 
a  detectable  size  because  the  Blake  threshold  exceeds  the  crevice  model 
threshold  in  the  near  saturation  region.  For  liquids  far  below 
saturation,  the  crevice  model  is  a  good  predictor  of  cavitation 
thresholds.  In  this  region  the  critical  radius  criterion  anoears  to  be 
a  necessary  and  sufficient  condition  not  only  for  bubble  p.ucleatlon  but 
also  for  hubble  detection. 

3-2.2  The  Effect  of  Surface  Tension 

The  functional  dependence  of  the  critical  radius  threshold  on  the 

liquid-gas  surface  tension  is  not  as  obvious  as  with  the  case  of  eas 

3  /+ 

pressure.  To  understand  this  dependence,  we  follow  the  lead  or  Crum  , 

who  used  an  empirical  relationship  discovered  by  Bargeman  and  Van  Voorst 
..  .  ^0 

Vader  ,  vi.z,  , 

c 

cos  a,.  *»  —  -  1 

E  V  (71) 

where  _  Is  the  equilibrium  contact  angle  and  c  is  a  constant.  For  many 
t 

nonpolar  solids  c  has  a  value  of  about  50  dyn/cm.  This  relationship  is 

valid  tor  surface  tensions  greater  than  aboiu  10  dyn/cm.  In  the  analysis 

given  by  Crum  as  well  as  the  one  to  follow,  it  is  assumed  that  Up=a^  and 

that  a  =a  -a  ,  where  a  is  the  hysteresis  contact  angle.  (Hysteresis  is 
RAH  M 
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what  gives  rise  to  the  asymmetrical  shape  of  droplets  sliding  down 

windows.)  Using  these  relationships,  the  crevice  model  becomes  a  model 

with  two  adjustable  parameters:  a  and  6.  In  the  analysis  of  this 

n 

model,  the  sum  of  these  two  parameters  +6)  was  held  constant.  The 

H 

value  of  the  sum  is  an  arbitrary  choice,  35°  giving  good  results. 

Figure  13  is  a  graph  of  threshold  versus  surface  tension.  Two  sets  of 

data  are  shown  corresponding  to  two  different  values  of  gas  pressure. 

34 

The  data  are  Crum's.  It  is  seen  that  the  agreement  between  theory  and 
experiment  is  quite  good.  It  should  be  pointed  out  that  the  receding 
contact  angle  criterion  offers  a  good  fit  as  well.  Thus,  while  these 
predictions  offer  no  direct  comparison  between  the  two  criteria,  they  do 
show  that  the  crevice  model  (in  either  form)  explains  the  dependence  on 
surfscc  Cession  quite  well* 


* 


t 
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3-2.3  The  Effect  of  Temperature 

Temperature  variations  affect  the  cavitation  threshold  through 
three  measureable  parameters:  gas  content,  surface  tension,  and  vapor 
pressure.  The  largest  variation  is  associated  with  the  temperature 
dependence  of  the  dissolved  gas  content;  however,  all  of  the 
dependencies  are  taken  into  account  in  the  following  analysis.  We 
examine  the  temperature  dependence  of  gas  pressure  first. 


Suppose  a  volume  of  water  is  in  equilibrium  with  a  gas  atmosphere 
of  pressure  P^0,  the  entire  system  being  at  temperature  T^.  Then 
according  to  Henry's  Law,  the  mole  fraction  of  the  gas  dissolved  in  the 


liquid,  X,  is  given  by  X°=Pg0/K°. 


K  is  called  Henry's  Law  constant  and 


t 


Graph  of  threshold  versos  surface  tension  .  The  data 
nd  the  solid  lines  are  calculated  from  equation  (31) 


62 


m 


s 


9 


9 


4 


4 


4 


is  a  function  of  temperature.  If  the  temperature  is  changed  to  T,  then 

gas  goes  into  or  out  of  solution  until  again  X**P  /K.  For  large  volumes 

g 

of  water,  the  amount  of  gas  which  diffuses  into  the  nucleus  is 

negligible.  Thus,  the  mole  fraction  can  be  considered  constant  and  the 

equilibrium  gas  pressure  at  a  temperature  T  becomes 

P  (T)  »  P  °  (72) 

g  g  K° 

where  the  temperature  dependence  is  indicated. 

In  order  to  evaluate  this  dependence ,  data  taken  from  the 

4  i 

International  Critical  Tables  for  the  variation  of  K  with  temperature, 
were  fit,  using  a  least  squares  method,  to  a  third  degree  polynomial, 
viz.  , 

R  -  (3-24xl0?)  +  (9.59xl05)  T  -  <2.17xl03)  T2  -  2.59  T3  .  (73) 

The  units  of  K  are  bar  per  mole  of  gas  per  mole  of  water.  T  is  measured 

in  degrees  C.  The  coefficient  of  correlation  of  this  fit  is  0.999. 

Figure  14a  is  a  graph  of  K  versus  temperature.  In  the  temperature  range 

of  interest,  the  fit  is  indeed  quite  good.  The  initial  value  of  K,  K° , 

is  calculated  from  a  knowledge  of  the  initial  gas  pressure  and  the 

temperature  at  which  the  sample  was  degassed. 

Although  analytical  expressions  for  the  temperature  dependence  cf 

surface  tension  exist,  in  this  analysis  we  chose  to  employ  a  least 

squares  fit  of  values  for  the  surface  tension  of  water  taken  from  the 

42 

CRC  Handbook  of  Chemistry  and  Physics  (56th  edition).  A  second  degree 
polynomial  least  squares  fit  was  used,  viz., 

Y(T)  «  75.65  -  0.142  T  -  (2.63xl0~4)  T2,  (74) 

Surface  tension  is  in  dyn/cm  and  temperature  in  degrees  C.  The 
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Figure  14a.  This  figure  illustrates  the  polynomial 
to  determine  Henry's  law  constant  as  a  function  of 
erature.  The  open  squares  are  data  from  reference  4 
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coefficient  of  correlation  is  0.999.  Figure  14b  displays  the  data  and 
least  squares  fit  for  the  surface  tension  dependence  on  temperature. 

The  final  temperatur  variation  to  be  considered  is  that  due  to 
vapor  pressure.  For  this  parameter,  an  analytical  expression  was  used. 
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The  expression  is  that  of  Osborne  and  Meyers  . 


B 


where 


log  P  -  A  +  -  + 

A  A 


CX  ,„DX 


10  -  1 


+  E 


10 


FXt 


(75) 


T  =  temperature  in  degrees  Celcius 

T.  =  absolute  temperature 
A 

X  =  T/  -  K 
A 

Y  =  374.11  -  T 
A  =  5.4267 


C  =  1.3869x10 
D  =  1.1965x10 
k  -  2.937xl05 
E  »  -0.0044 
F  »  -5.715x10 


-4 


-11 


-3 


B  =  -2005.1 

According  to  Eisenberg  and  Rauzmann^3,  this  expression  is  accurate  to 

4 

within  five  parts  in  10  . 

The  cavitation  threshold  is  plotted  as  a  function  of  temperature 
for  a  range  of  dissolved  gas  pressures  in  figure  15.  The  open  squares 
represent  Gateway's  data  and  the  solid  squares  represent  Crum's  data. 
The  solid  lines  represent  the  threshold  as  predicted  by  the  critical 
radius  criterion  for  6=12.85°  and  8+a=35°.  The  data  have  been 

n 

normalized  to  fit  the  theory  at  room  temperature  (23°  C).  The  variation 
with  temperature  is  predicted  well  by  the  crevice  model  employing  the 
critical  radius  criterion.  Although  there  are  no  data  to  test  the 

o 

theory  above  50  C,  the  threshold  must  approach  the  vapor  pressure  at 
the  boiling  point.  This  behavior  is  suggested  by  the  "tail'1  at  high 
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Figure  15.  Graph  of  threshold  versus  temperature.  The  data  are 
Galloway's  (Q)  and  Crum's  (■)  and  have  been  normalized  to  fit  thi 
theory  at  room  temperature.  The  solid  lines  are  calculated  with 
equation  (31)  using  6=12.85°  and  a. =22, 15°. 
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temperatures. 

3-2.4  Comments  on  Sensitivity 

It  was  shown  that  the  critical  radius  approach  to  the  crevice  model 
provides  a  good  explanation  of  acoustic  cavitation  thresholds.  However, 
all  of  these  predictions  have  beer  made  a_  posteriori.  That  is,  values 
of  0  and  $  were  chosen  to  provide  the  best  fit  to  the  data.  For  the 
majority  of  the  analysis,  8  was  chosen  to  be  12.85°  and  <j>  “35°.  The 
exception  was  the  prediction  of  threshold  versus  surface  tension  (figure 
13)  for  P^=0.19  bar.  For  this  line  8-10.8°,  yet  $  remained  equal  to 
35°.  The  model  is  very  sensitive  to  changes  in  B*  Figure  16a  shows 
theoretical  predictions  of  threshold  versus  surface  tension.  In  each 
line,  the  parameter  d  is  constant  and  c^ual  to  35°.  For  the  middle 
line,  0^=22.15°  and  8-12.85°.  These  are  the  same  parameters  for  the 
P^»0.47  bar  prediction  in  figure  13.  The  upper  curve  is  for  a^-21.150, 
8-13.85°;  the  lower  curve  is  for  a  -23.15°,  8  -11.85°.  Figure  16b  is 

n 

similar  to  figure  16a,  however,  in  this  figure  a  -22.15°  for  all  the 

H 

curves  and  8-12.85  +  1°  for  Che  middle,  upper  and  lower  curves 
respectively.  In  figure  16c,  8-12.85°  and  o.  ^-22.15  +_  1°  for  the  middle, 
upper  and  lower  curves ,  respectively. 

It  is  seen  that  varying  a  has  little  effect.  However,  variation 
oi  8  produces  a  large  change  in  predicted  thresholds.  This  is 
unfortunate.  It  means  that  extremely  accurate  a^  priori  predictions  of 
cavitation  thresholds  are  n>  :  possible.  However,  the  purpose  of  this 
research  is  not  to  predict  cavitation  thresholds,  per  se ,  but  to 
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Figure  16a.  Graph  of  threshold  versus  surface  tension 
illustrating  the  sensitivity  of  equation  (31)to  variou 
parameters . 


Figure  16c.  Graph  of  threshold  versus  surface  tension 
illustrating  the  sensitivity  of  equation  (31)  to  changes 
in  various  parameters. 


71 


identify  possible  candidates  for  the  cavitation  nucleus.  The  degree  to 
which  the  a  posteriori  predictions  agree  with  measured  values  indicate 
that  the  crevice  model  is,  indeed ,  a  plausible  r.ucleation  model  for 
acoustic  cavitation.  In  addition,  even  in  the  wors t  case ,  using  typical 
values  of  the  parameters,  a_  priori  predictions  of  the  cavitation 
threshold  can  be  made  that  disagree  with  measured  thiesholds  by  no  more 
than  a  factor  of  two- 

3-3  Application  cf  the  Crevice  Model  to  Diffusion  Cavitation 

The  analysis  presented  in  this  section  is  unique;  ic  is  the  first 

attempt,  that  we  know  of,  at  applying  thj'  crevice  model  to  diffusion 

cavitation.  Reference  will  be  made  to  rhe  diffusion  cavitation  data 
26 

published  by  Yount  ,  an  example  of  which  was  reproduced  as  figure  4, 
Yount  has  explained  his  data  in  terms  of  the  varying-peimeability 
orga.  kin  model.  The  variation  of  the  permeability  is  necessary  to 

•  c  change  in  slope  of  the  isopleths  for  the  various  numbers  of 

bv.  ' ibis  change  in  slope  is  something  which  Yount  claims  is 

inexplicable  in  terms  of  the  crevice  model.  The  version  of  the  crevice 

model  to  which  Yount  is  referring  to  is  the  previous  formulation.  An 
explanation  of  this  phenomenon  can  be  provided  employing  the  revised 
crevice  model. 

* 

As  pointed  out  in  the  previous  chapter,  the  concept  of  P  crus^  is 
analogous  to  that  of  a  critical  size  of  a  crevice.  Yount  also  ur.es 

X 

P  cruc.^  t0  mean  the  crushing  pressure  at  which  the  organic  skin  becomes 

•k 

impermeable.  Drawing  an  analogy  between  the  two  definitions  of  P  ,  , 
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we  can  say  that  the  crevice  model  analogy  of  a  permeable /impermeable 

transition  is  the  interface  reaching  the  advancing  contact  angle. 

* 

Equation  f 63 )  was  derived  for  crushing  pressures  greater  than  P  ,  ; 

crush 

i.e.,  it  is  analogous  to  a  process  in  Yount's  impermeable  region.  For 

* 

crushing  pressures  less  than  P  ,  ,  equation  (63)  becomes 

crush 


P 


cos(aR-g) j 


ss  crush  |cos(ct  -0) 


(1 +  nRtan3) 3  . 


(76) 


f 
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The  only  difference  between  equations  (63)  and  (7b)  is  that  in  equation 
(76),  the  interface  is  at  an  angle  a  and  not  the  advancing  contact  angle 
a  .  Moreover,  9G°  <  a  -  8  <_  a-S  ,  which  means  that  |cos(a-8  )|  <_ 

|cos(aA  -d  )j.  Hence,  d(Pfi6)/d(Pc  sh)  8reater  when  the  crushing 

if 

pressure  is  less  Chan  P  ..  This  result  is  consistent  with  Yount's 

ccusn 

findings. 

fn  the  experiments  described  by  Yount,  bubbles  were  formed  in 
gelatin  samples  which  had  an  average  surface  tension  of  about  53  dyn/cm. 


Using  equation  (71),  one  finds  that  the  advancing  contact  angle  for  this 

surface  tension  is  around  93°.  It  can  be  seen  from  geometrical 

v 

considerations  -hat  for  stability  a  >y+B>  a  .  Therefore  if 

A  '■  K 

a  “  93°,  the  crevice  half~angle  must  be  between  zero  and  three  deg..-. 

A 

Crevices  with  half-angles  greater  than  3°  are  unstable  and  become 

deactivated  (thoroughly  wetted).  For  the  calculations  to  follow,  p  was 

assumed  to  be  1°;  the  crevice  has  become  a  narrow  crack.  For  such  a 

narrow  crevice  with  an  advancing  contact  angle  close  to  90°,  one  would 

expect  that  the  hysteresis  angle  would  be  small;  assume  it  to  be  3J.  A 
*  26 

typical  value  of  P  .is  about  8  bars.  Solving  equation  (47)  for 
f  crush 
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the  product  hA  and  using  these  values  for  the  various  parameters,  one 
finds  that  hA  is  about  6  nm.  One  possible  interpretation  of  this  result 
is  that  the  nucleus  is  located  at  the  bottom  of  a  deep,  narrow  crevice. 

A 


For 
d(P  )/d(P 


F  <  P 

crush  crush' 


one  finds 
between  0,5  and 


from  equation  (76)  that 

A 


P  >  P 

crush  —  crush' 


..  l  ,)  is  between  0,5  and  1;  for 
8  S  C lU 8  n 

* 

d(P„)/d(F  crufl|l>  about  0*5.  These  values  are  reasonably  consistent 
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with  measured  values. 
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The  main  difference  is  that  the  data  show  a 


range  of  3lope3  for  P  .  >  P  ,  from  about  0.2  to  1.  The  crevice 

crush  crush 

model  result  is  within  this  range.  This  discussion  is  illustrated  in 
figure  17.  The  predictions  shown  in  this  figure  are  reminiscent  of 
those  in  figure  4. 

The  results  of  these  calculations  are  very  sensitive  to  the  chosen 

parameters  and  more  work  needs  to  be  done  on  these  concepts.  The  size 

of  the  nucleus  involved  with  these  processes  is  very  small  and  it  is 

possible  that  some  phenomenon  not  taken  into  account  in  this  model  may 

well  become  important  at  these  small  sizes.  However,  two  points  have 

been  made.  First,  the  crevice  model  can  explain  the  change  in  slope  of 

the  isopleths.  Second,  the  calculated  values  of  the  slopes  are 

reasonable.  Perhaps  one  final  comment  it  appropriate.  Yount  and  others 

have  carried  out  an  extensive  microscopic  investigation  of  bubble 
44 

formation  nuclei,  about  which  more  will  be  discussed  later.  The 
investigators  make  the  statement  "In  no  case  have  we  observed  a  gas 
phase  embedded  in  a  crevice,  as  would  be  required  by  the  crevice 
model..."  This  absence  is  exactly  what  one  would  expect  according  to 


( 


the  argument  given  above,  which  showed  that  the  gas  phase  should  be  very 
small  and  located  at  the  bottom  of  a  deep  narrow  crevice. 

As  has  just  been  shown,  the  critical  radius  approach  to  the  crevice 
model  offers  an  explanation  of  results  of  both  acoustic  and  diffusion 
cavitation  experiments.  There  are,  however,  other  nucleation  models  and 
they  are  taken  up  next. 


Chapter  4 


Other  Candidates  for  the 
Cavitation  Nucleus 

4-1  Introduction 

The  previous  chapters  have  been  devoted  to  the  critical  radius 

approach  to  the  crevice  model.  It  was  shown  that,  based  upon  its  a 

posteriori  predictive  powers,  it  is  a  plausible  nucleation  model.  There 

are,  however,  alternative  models  and  this  chapter  is  devoted  to  them. 

It  was  stated  in  the  introduction  that  only  two  other  nucleation  models 
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are  cor«idered  feasible.  One  is  the  ionic  skin  model  of  Akulichev  and 
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the  other  is  the  organic  skin  model  most  recently  advanced  by  Yount. 

In  this  chapter,  these  models  will  be  analyzed  to  determine  whether  or 
not  they  really  are  plausible  models. 

4-2  Ionic  Skin  Model 

The  ionic  skin  model  proposed  by  Akulichev  was  reviewed  in  the 
first  chapter.  That  review  is  adequate  for  this  analysis  and  the  reader 
should  refer  back  to  it  if  necessary.  The  strongest  evidence  that 
Akulichev  presents  in  support  of  his  model  is  the  increase  in  threshold 
with  decreasing  ion  concentration.  This  behavior  must  exist  if  his 
model  is  feasible.  It  is  to  this  behavior  that  we  direct  our  attention. 

Two  points  become  evident  upon  investigating  Akulichev' s  work. 
First,  the  lower  limit  of  the  concentrations  used  in  his  experiments  was 


about.  1  ramol/'l.  In  fact,  he  published  only  two  data  points  below  this 
concentration,  one  at  about  0.6  mmol/1  and  the  other  about  0.3  inmol/1. 
Second,  his  conclusions  about  thresholds  for  concentrations  below 
1  mmol/1  are  based  upon  extrapolations  of  trends  present  at  higher 
concentrations,  not  upon  experimental  measurements.  (Refer  to  figure  1.) 
In  an  effort  to  extend  this  lower  limit  and  experimentally  determine  the 
influence  of  the  dissolved  ion  concentration  on  the  cavit  .tion 
threshold,  we  performed  a  series  of  measurements  of  the  effect  of  ions 
on  the  cavitation  threshold.  The  details  of  this  experiment  are 
presented  in  the  appendix;  only  those  details  necessary  for  the  present 
discussion  will  be  mentioned  here. 

Figures  18a  and  18b  represent  typical  results  of  the  experiment. 
It  should  be  pointed  out  that  Kero  on  the  abscissia  corresponds  to  a 
concentration  of  1  mmol/1.  Tne  data  in  the  figures  have  been  normalized 
to  the  threshold  value  at  the  lowest  concentration.  Tne  experimental 
scatter  is  such  that  no  strong  conclusions  can  be  drawn  about  the 
possible  fine  structure  present  iri  the  data.  however,  one  conclusion 
can  be  made.  The  threshold  does  not.  increase  racnotonically  with 
decreasing  concentration.  On  the  contrary,  in  every  set  of  measurements 
made,  the  trend  is  for  the  threshold  to  decrease.  It  is  true  that  no 
direct  comparison  of  these  data  and  Akulichev’s  is  possible  because  the 
experiments  were  performed  at  different  dissolved  gas  concentrations. 
However,  it  is  not  at  all  clear  why  this  difference  in  gas  content 
should  alter  the  results  so  as  to  produce  the  opposite  behavior.  It 
may,  therefore,  be  concluded  that  the  nuclei  of  the  cavitation  detected 


Figure  18a.  Graph  of  normalized  cavitation  threshold  versus  the 
log  of  the  dissolved  ion  concentration  for  KI . 


normalized  cavitation  threshold  versus  the 
concentration  for  Nal. 


in  our  experiments  were  not  ionically  stabilized. 

Although  not  statistically  significant,  there  is  a  recurring  trend 
in  all  of  the  data  for  the  threshold  to  reach  a  maximum  value  at  a 
concentration  of  about  1  mmol/1.  This  behavior  is  reminiscent  of  the 
trend  in  Akullchev'o  data.  It  is  an  interesting  observation  and  it 
makes  theoretical  analysis  more  difficult.  This  trend  leads  to  a  second 
conclusion  that  can  be  drawn  from  the  ion  data.  The  dissolved  ion 
concentration  does,  in  fact,  influence  the  cavitation  threshold.  This 
influence  is,  perhaps,  not  a  simple  one  and  it  appears  that  it  is  not 
attributable  to  a  purely  ionic  skin.  The  trend  in  the  data  may  indicate 
a  transition  region,  above  which  one  type  of  nucleus  is  predominant  and 
below  which  a  second,  different,  type  of  nucleus  is  the  predominant 
cause  of  the  cavitation.  For  example,  the  cavitation  detected  at 
concentrations  above  1  mmol/1  may,  in  fact,  be  due  to  Akulichev's 
nuclei.  Decreasing  the  concentration  would  increase  the  threshold  of 
these  nuclei.  However,  the  decrease  in  concentration  may  decrease  the 
threshold  of  some  other  type  of  nucleus  so  that  it  is  able  to  produce  a 
bubble  at  a  lower  threshold  than  an  Akulichev  nucleus.  This  statement 
is  speculative  but  plausible. 

4-3  Varying-Permeability  Model 

In  section  3-3,  the  crevice  model  was  applied  to  diffusion 

cavitation,  a  process  which  had  been  explaineu,  almost  exclusively  in 

2  8 

terms  of  the  organic  skin  model  .  In  this  section,  the  organic  skin 
model  will  be  applied  to  acoustic  cavitation,  a  process  which  has  been 
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explained  well  in  terras  of  the  crevice  model.  The  dependence  on 
dissolved  gas  pressure  will  be  discussed  first. 

In  order  to  determine  the  variation  in  the  acoustic  cavitation 
threshold  with  dissolved  gas  pressure  based  upon  an  organic  skin 
nucleus,  it  must  first  be  determined  how  such  a  nucleus  is  affected  by 


the  degassing  process.  The  degassing  of  a  liquid  is  generally 
accomplished  in  the  following  manner.  The  atmosphere  surrounding  the 
liquid  is  evacuated  to  a  certain  pressure,  say  P^  ("d"  for  degassing 
pressure).  This  is  equivalent  to  the  decompression  phase  of  Yount’s 
experiments.  The  liquid  comes  to  equilibrium  with  the  atmosphere  so 
that  the  dissolved  gas  pressure  equals  P^.  At  this  point  the  ambient 
pressure  is  returned  to  atmospheric  pressure,  which  is  equivalent  to  the 
compression  phase  of  Yount's  experiments.  The  measurements  are  then 
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made  before  the  liquid  has  enough  time  to  come  back  into  equilibrium 
with  the  atmosphere.  The  process  just  described  is  a  Yount  experiment 


in  reverse;  the  decompression  phase  occurs  first,  then  the  equilibratioa 
phase,  and  finally  the  compression  phase.  Therefore,  in  order  to 
determine  the  effect  of  gas  pressure,  the  effect  that  the  degassing 
process  has  on  the  radius  of  the  nucleus  must  be  determined.  Since  the 
changes  in  pressure  are  always  less  than  one  atmosphere,  the  skin  can 
always  be  considered  permeable. 

Assume  that  a  varying-permeability  (V-P)  nucleus  with  an  initial 
radius  Rq  is  degassed  to  a  pressure  .  This  decompression  is  governed 
by  equation  (19c)  of  reference  26,  viz.. 
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(The  notation  has  been  altered.)  is  called  the  skin  compression. 

The  decompression  results  in  the  liquid  being  supersaturated,  so  the 


nucleus  will  grow  bv  diffusion  if  the  final  radius  R,  exceeds  R  . 

d  eg 

Using  this  criterion  along  with  the  preceeding  equation,  one  finds  that 
the  largest  radius  which  a  nucleus  can  have  and  survive  degassing  (i.e. 


will  not  grow  by  diffusion)  is 
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According  to  the  V-P  model,  the  following  relationship  exists  among  Y  , 


Y  and  R  : 
o 


Yc  =  Y  +— 


(The  variable  b  is  the  same  as  beta  in  reference  26.)  b  is  a  variable 

which  depends  upon  thermodynamic  considerations  but  is  constant  for  a 

given  sample.  For  the  experiments  discussed  by  Yount  it  has  a  value  of 
5  2 

about  1.43  x  10  N/m  .  Substituting  equation  (79)  into  (78)  gives 


R  =  2  Y 
o 
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Following  the  decompression,  the  liquid  eventually  comes  to 
equilibrium  with  the  reduced  atmosphere.  It  is  not  clear  exactly  what 
the  response  of  the  nucleus  will  be  to  this  equilibration.  However, 
two  possible  extreme  cases  can  be  discussed.  The  first  case  is  one  in 
which  the  equilibrium  radius  ,  R^,  is  the  same  as  the  degassed  radius, 
R^ .  In  order  for  this  to  be  the  actual  case,  enough  additional 
surfactant  molecules  must  be  adsorbed  onto  the  skin  t:  stabilize  the 
skin  at  this  radius.  The  second  extreme  case  is  one  in  which  the 


nucleus  shrinks  back  to  its  initial  radius,  R  .  In  this  case  it  is 

o 
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assumed  that  no  additional  surfactant  molecules  attach  themselves  to 
the  skin.  Moreover,  one  would  not  expect  any  molecules  to  leave  the 
skin  during  the  decompression.  Hence,  since  the  skin  consists  of  the 
same  number  of  molecules  the  nucleus  should  be  stabilized  at  the  same 
radius,  namely,  R  . 

In  either  case,  during  the  subsequent  compression  back  to  P  ,  the 

26 

nucleus  should  obey  the  relationship 
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where  is  the  final  radius  after  compression,  is  the  equilibrium 

radius  (either  R  or  R , )  ,  ard  P  is  the  dissolved  gas  pressure  in  the 
o  d  g 

liquid  and  equal  to  P^.  Comparing  equations  (77)  and  (81),  it  is  easily 

seen  that  if  the  first  case  Is  true  (R  =*  R, )  then  R  =  R  .  If  the 

e  d  r  o 

second  case  is  true  then  it  can  be  shown  from  equations  (79)  and  (81) 


t.  =  R  1 
t  o 


P  -  P 

O  £ 


The  compressed  nucleus  will  be  supersaturated  so  gas  will  diffuse  out  of 

it;  however,  this  diffusion  does  not  result  in  a  change  in  radius  (see 

equation  '6  of  reference  26).  Now  that  the  response  of  a  V— P  nucleus  to 

the  degassing  process  is  known,  it  roust  determined  how  such  a  nucleus 

will  react  to  an  applied  acoustic  field.  Previously,  it  was  stated  that 

a  V~P  nucleus  will  grow  during  degassing  if  its  radius  exceeds  R  . 

This  criterion  is  consistent  with  the  derivation  of  the  V-P  model  (see 

equation  17a  of  reference  26).  There! ore,  one  would  expect  that  such  a 

nucleus  will  grow  in  an  acoustic  field  when  the  radius  exceeds  R  (if 
°  c  v 


P,  <  P  -P  ).  This  criterion  is  given  by  the  Blake  threshold  (see 
A  o  v  1 

equation  (69)  of  Chapter  3).  Figure  19a  shows  the  Blake  threshold  as  a 
function  of  dissolved  gas  pressure  for  the  two  extreme  cases  discussed 
above.  It  can  be  seen  that  the  threshold  is  relatively  insensitive  to 
variations  in  gas  pressure  and  that  the  difference  between  the  two  cases 
is  negligible.  However,  there  is  a  hint  that  the  proper  trend  is 
present. 

There  is  one  large  caveat  in  this  discussion;  it  is  the  choice  of 
the  value  of  b.  For  the  results  shown  in  figure  19a,  the  value  of  b  was 
that  computed  in  reference  26.  That  is,  it  was  based  on  measurements 
made  on  gelatin  samples  having  a  surface  tension  of  about  51  dyn/cm. 
However,  the  data  presented  in  figure  \9a  were  taken  in  water  having  a 
surface  tension  of  about  79  dyn/cm.  The  dependence  of  b  on  surface 
tension  is  given  in  equation  (39b)  of  reference  26, 

b  »  2.80  x  10^Y  dyn/cra.  (83) 
For  Y  “  70  dyn/cm,  b  «■  1.96  bar.  Figure  19b  is  a  graph  of  threshold 
versus  gas  pressure  for  b  »  1.96  bar.  Comparing  figures  19a  and  19b,  it 
is  seen  that  this  modification  snakes  the  prediction  worse  instead  of 
better. 

Upon  examination  of  equation  (82)  it  is  seer  that  if  b  is 
increased,  R^  increases  as  well;  a  larger  R^  results  in  a  lower 
threshold.  Therefore,  increasing  b  from  1.43  to  1.96  bar  results  in  the 
threshold  decreasing.  In  order  for  the  V-P  model  to  predict  the 
acoustic  cavitation  threshold  as  a  function  of  cissolved  gas  pressure,  b 


Figure  19a.  Graph  of  the  cavitation  threshold  versus  dissolved 
gas  pressure  for  a  V-P  nucleus  with  b=1.43  bar.  (Refer  to  j.ig- 
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nucleus  with  b=L.OO  bar.  (Refer  to  fig- 
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must  be  less  than  1.43  bar,  not  greater.  The  lower  limit  on  b  is  1  bar. 
If  b  were  less  than  1  bar,  equation  (80)  would  predict  negative  radii 
for  seme  values  of  Pgg.  Figure  19c  represents  the  threshold  dependence 
on  gas  pressure  for  b  *  1  bar.  The  agreement  is  much  better,  however, 
the  means  by  which  the  agreement  was  improved  is  not  consistent  with  the 
V-P  model. 

From  this  discussion,  either  of  two  conclusions  may  be  drawn.  Hie 
first,  and  most  severe,  conclusion  that  may  be  drawn  is  that  the  V-P 
model  is  not  adequate.  A  second,  and  perhaps  more  realistic,  conclusion 
is  that  the  experimental  confirmation  of  the  V-P  model  has  been 
inadequate.  The  only  measurements  to  which  the  V-P  model  has  been 
applied  have  involved  diffusion  cavitation  occurring  in  gelatin.  The 
numerical  factor  (2.80)  present  in  equation  (83)  is  based  upon  a  set  of 
measurements  made  on  a  single  sample  of  gelatin.  It  is  possible  that 
this  factor  varies  with  the  experimental  conditions  and  is  not  a 
"universal"  constant.  This  possibility  could  account  for  the 
discrepancy  present  in  figure  19. 

Aside  from  the  limited  experimental  confirmation  of  the  V-P  model, 
its  complexity  leaves  it  cumbersome  to  apply.  For  instance,  consider 
the  effect  that  a  variation  in  surface  tension  might  have  on  the 
nucleus.  Since  the  model  depends  upon  the  chemical  potential  of  the 
surfactants,  It  is  not  at  all  clear  that  if  different  types  of 
surfactants  are  used  to  alter  the  surface  tension,  the  threshold  would 
not  be  affected  in  different  wv. . s.  Even  the  model's  chief  advocate  is 
not  certain  how  the  addition  of  surfactants  affect  the  nucleus.  In 
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reference  28,  the  statement  is  made  "Another  standard  method  [of 

increasing  the  threshold  fcr  bubble  formation]  is  to  add  detergents  or 

other  substances  which  presumably  operate  by  changing  [bl".  The  effect 

of  temperature  variation  on  a  V-P  nucleus  is  an  equally  complex  matter. 

In  fact,  not  even  a  qualitative  description  of  the  temperature 

dependence  is  offered  in  reference  28,  only  the  statement  "A  dependence 

on  temperature  is  also  predicted". 

The  discussion  in  this  section  may  be  summarized  as  follows. 

Although  the  vsrying-perraeability  model  is  elegant  and  explains 

diffusion  cavitation  in  gelatin,  its  application  to  acoustic  cavitation 

is  cumbersome  at  best.  Since  these  nuclei  have  been  observed  using 

42 

microscopic  techniques  ,  their  existence  is  difficult  to  deny. 
However,  the  model  needs  to  be  extended  to  explain  acoustic  cavitation 
processes. 
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Summary,  Conclusions  and 
Topics  for  Further  Study 


5-1  Summary  of  the  Dissertation 

As  stated  at  the  outset,  the  purpose  of  the  research  discussed  in 


this  dissertation  was  to 

advance 

the  state  of 

knowledge 

of 

Che 

cavitation  nucleus.  The 

advances 

made  toward 

this  goal 

may 

be 

summarized  as  follows; 

1)  The  crevice  model  was  revised  in  terms  of  a  more  physically 
realistic  nuclcsticn  criterion.  This  new  criterion  was  shown  rn  be  both 
a  necessary  and  sufficient  condition  for  cavitation  nucleation. 

2)  The  predictions  of  the  acoustic  cavitation  threshold  based  on 
this  revised  model  (.called  the  critical  radius  approach  to  the  crevice 
model)  were  found  to  be  in  good  agreement  with  experimental  results. 

3)  The  revised  crevice  model  was  applied  to  diffusion  cavitation 
(the  first  such  application)  and  it  was  shown  to  predict  results 
previously  explained  only  in  terms  of  the  varying-permeability  model. 

4)  It  was  shown  that  in  liquids  with  low  surface  tensions, 
nanometer-sized  nue’ei  are  capable  of  producing  cavitation. 

5)  The  first  extensive  measurements  of  the  acoustic  cavitation 
threshold  as  a  function  of  dissolved  ion  concentration  were  made  for 
concentrations  less  than  1  mtaol/1. 

6)  These  measurements  indicate  that  the  cavitation  observed  could 
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not  have  been  produced  entirely  by  a  purely  ionic  skin  as  proposed  by 
Akulichev. 

7)  The  varying  permeability  model  was  applied  to  acoustic 
cavitation  (another  unique  application).  It  was  found  that  although  the 
model  could  explain  the  increase  in  threshold  with  decreasing  gas 
content,  certain  inconsistencies  arose  which  indicated  that  the  model 
should  be  developed  further  to  facilitate  its  application  to  acoustic 
cavitation. 

5-2  Conclusions 

In  this  dissertation,  three  different  cavitation  nucleation 
mechanisms  were  considered.  Based  upon  the  findings  of  this  work,  v;e 
conclude  that  the  critical  radius  approach  to  the  crevice  model  is  the 
most  satisfactory  of  these  models.  Using  this  model,  the  threshold  for 
acoustic  cavitation  can  be  accurately  predicted  as  a  function  of  various 
parameters.  In  addition,  some  aspects  of  diffusion  cavitation  are  seen 
to  be  natural  consequences  of  the  crevice  model.  We  should  stress, 
however,  that  this  conclusion  is  based  on  measurements  made  in  water. 
It  should  not  be  assumed  that  these  findings  can  be  extrapolated  to 
other  fluids  without  some  detailed  analysis. 

Cavitation  is  a  vastly  complicated  phenomenon  and  to  conclude  that 
a  single  type  of  nucleus  accounts  for  all  cavitation  processes  is  almost 
certainly  in  error.  We  do  feel  that  the  critical  radius  approach  to  the 
crevice  model  is  the  most  satisfactory  theory  at  present.  Indeed,  the 
crevice  model  may  account  for  the  largest  part  of  the  entity  referred  to 


t 


as  "the"  cavitation  nucleus.  However,  the  fact  remains  that  there  is  no 

theoretical  reason,  nor  any  conclusive  experimental  results,  that 

discounts  the  existence  of  the  other  two  proposed  nuclei.  (In  fact, 
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organic  skin  nuclei  have  been  ohserved  by  Yount.  )  These  models 
should  be  developed  further.  As  was  stated  earlier,  we  have  advanced 
the  state  of  knowledge  concerning  the  nucleation  of  cavitation  in  water. 
The  next  step  to  be  taken  should  be  toward  uniting  these  models  into  "a" 
cavitation  nucleation  model. 

5-3  Topics  for  Further  Study 

As  the  research  presented  in  this  dissertation  progressed,  it 
became  clear  that  some  aspects  of  this  research  deserved  further  study. 
In  the  previous  section  some  of  these  topics  were  mentioned .  In  this 
section  we  will  mention  some  others. 

1)  Here  measurements  need  to  be  taken  in  different  fluids.  In 
particular,  biological  fluids  are  presently  in  need  of  extensive 
investigation  because  of  the  possible  hazards  associated  with  the  use  of 
ultrasound  in  medicine. 

21  The  results  of  the  measurements  of  threshold  as  a  function  of 
ion  concentration  are  not  easily  explained  in  terms  of  our  present 
knowledge.  The  explanation  of  these  results  has  the  potential  of 
advancing  our  knowledge  of  nucleation  processes  significantly. 

3)  The  modification  of  the  experimental  apparatus  used  to  make 
the  ion  measurements  ( cc  be  discussed  in  the  appendix)  should  make  it 
possible  to  verify  the  existence  of  the  surfactant  or  ionic  skin  nuclei. 
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With  these  modifications  it  should  be  possible  to  remove  the  vast 
majority  of  the  existing  cavitation  nuclei.  Bubbles  could  then  be 
introduced  into  the  liquid  (perhaps  through  electrolysis)  in  the 
presence  of  surfactants  or  ions.  Since  some  of  these  bubbles  should 
stabilize,  the  result  should  be  a  (significantly)  lower  caviation 
threshold . 

4)  A  further  modification  of  the  apparatus  would  make  it  possible 
to  perform  threshold  measurements  with  ambient  liquid  pressures  close  to 
the  vapor  pressure.  To  the  author's  knowledge,  this  type  of  measurement 
has  never  been  performed.  The  results  of  such  measurements  are  easily 
predictable  in  terms  of  the  revised  crevice  model,  and  therefore  this 
measurement  should  be  a  good  test  of  the  model. 

5)  Another  topic  for  further  study  is  the  measurement  of  the 
cavitation  threshold  as  a  function  of  dissolved  gas  pressure  for  ambient 
pressures  greater  than  atmospheric.  Such  an  experiment  may  prove  useful 
in  understanding  more  fully  the  reason  for  the  discrepancy  between  the 
predicted  threshold  based  on  equation  (3) )  and  the  measured  threshold 
for  gas  pressures  near  saturation  (refer  to  figure  11a). 

6)  The  last  topic  we  mention  is  of  a  different  sort.  Assume  that 
molecular  skin  nuclei  are  confirmed  to  exist.  Certain  types  of 
molecules  should  be  better  stabilizers  than  others.  Nuclei  composed  of 
good  stabilizers  may  be  relatively  rugged  and  able  to  exist  in  harsh 
liquid  conditions.  These  nuclei  would  be  long  lived"seeds"  which  could 
be  activated  at  a  certain  time  for  a  certain  purpose.  The  development 
of  such  technology  may  prove  to  be  quite  useful. 
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Append ix 


Experimental  Determination  of  the  Variation  of  the  Acoustic. 
Cavitation  Threshold  as  a  Function  of  Dissolved  Ion  Concentration 

A-l  Introduction 

In  this  appendix  the  experimental  aspects  of  the  measurement  of  the 
acoustic  cavitation  threshold  of  water  as  a  function  of  the  dissolved 
ion  concentration  of  the  water  will  be  discussed.  This  discussion  will 
be  presented  in  four  parts:  the  purpose  of  the  experiment,  the 
experimental  apparatus,  the  experimental  procedure,  and  the  results  of 
the  experiment.  The  conclusions  drawn  from  the  expe  iment  were 
presented  in  section  4-2. 

A-2  Purpose  of  the  Experiment 

One  purpose  of  the  experiment  was  to  attempt  to  verify  the  results 

20 

of  Akulichev's  experiment  since  his  results  are  the  only  strong 
evidence  for  the  existence  of  an  j.onically  stabilized  nucleus.  In 
addition,  our  apparatus  made  it  possible  to  extend  the  range  of  his 
measurements  to  lower  ion  concentrations.  This  range  of  concentrations 
(below  1  mraol/1)  had  not  previously  been  investigated,  and  yet,  the 
behavior  of  the  threshold  at  low  concentrations  is  a  strong  test  of  the 
ionic  skin  hypothesis. 

A-3  Experimental  Apparatus 

The  apparatus  used  in  making  these  measurements  was  developed  over 
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a  period  of  several  years  by  our  research  group.  It  is  designed  to  be 

versatile  so  that  it  can  be  used  to  make  threshold  measurements  for  a 

wide  range  of  fluid  conditions.  Since  the  system  is  described  in  some 
45 

detail  elsewhere,  only  a  brief  description  will  be  given  here. 

The  apparatus  is  composed  of  two  major  groups:  the  sample 

preparation  system  and  the  data  acquisition  system.  The  sample 
preparation  system  is  shown  schematically  in  figure  20a.  A  typical 
sample  preparation  takes  place  as  follows.  Prefilttred  distilled  water 
is  introduced  into  the  reservoir.  The  water  is  pumped  from  the 
reservoir  and  may  be  passed  through  either  the  organic  removal  column  or 
the  ionic  removal  column,  or  both.  The  organic  removal  column  is  a 
Barnstead  Organic  Removal  Cartridge  (part  number  D8904)  and  lowers  the 
concentration  of  long  chain  organic  molecules  to  less  than  one  part  per 
million.  The  ionic  removal  column  is  a  Barnstead  Ultrapure  Cartridge 
(part  number  D8902)  and  produces  water  with  a  conductivity  of  less  than 
one  micromho  and  a  pH  of  between  6.8  and  7.2.  The  specific  conductance 
of  the  water  is  measured  with  a  Cole  Parmer  conductivity  meter  (model 
number  R-1491-20).  The  resolution  of  the  meter  is  0.1  micromho  and  the 
accuracy  is  0.2  micromho.  The  water  then  passes  through  an  in-line 
teflon  filter  to  remove  solid  contaminants.  The  filter  size  can  be 
chosen  to  suit  the  requirements  of  the  experiment.  After  the  filter, 
tiie  water  circulates  through  a  spherical  quartz  flask  in  which  the 
acoustic  field  is  produced,  and  then  back  to  the  reservoir.  Except  when 
passing  through  the  organic  and  ionic  removal  columns,  the  water  is  in 
contact  with  either  teflon,  glass  or  quartz.  Tn  this  way,  the 
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possibility  of  contamination  of  the  sample  is  reduced.  The  capacity  of 
the  system  is  about  two  liters. 

The  gas  content  of  the  system  is  controlled  in  the  following 
manner.  The  glass  reservoir  is  evacuated  to  a  specified  pressure.  A 
gas  (usually  argon)  i3  bubbled  through  the  liquid  in  the  reservoir  as 
the  liquid  is  circulated.  The  agitation,  so  produced,  aids  in  both 
removing  the  dissolved  air  and  saturating  the  liquid  with  the  gas  at  the 
specified  pressure.  Before  entering  the  reservoir,  the  gas  is  passed 
through  a  filter  to  prevent  the  introduction  of  solid  contaminants. 

The  quartz  flask  is  enclosed  in  a  light  tight  box  which  contains 
two  large  liquid  heat  sinks.  The  heat  sinks  are  maintained  at  a 
constant  temperature  to  aid  in  controlling  the  temperature  of  the  sample 
during  the  data  acquisition  phase.  Also,  the  reservoir  is  surrounded  by 
a  cooling  jacket  which  is  used  to  adjust  the  temperature  of  the  sample. 
The  temperatures  of  the  sample  and  the  heat  sinks  are  monitored  with 
thermistors. 

The  data  acquisition  system  is  depicted  schematically  in  figure 
20b.  The  acoustic  field  is  generated  by  two  cylindrical  PZT  transducers 
cemented  onto  the  sides  of  the  quartz  flask.  They  are  driven  by  the 
amplified  signal  of  the  oscillator  In  a  resonant  mode  near  60  kHz.  The 
resonance  frequency  is  determined  by  monitoring  the  output  of  a  pill 
tr?  isducer  (which  is  also  cemented  on  the  side  of  the  flask)  on  an 
oscilloscope.  The  amplitude  of  the  driving  signal  is  modulated  by  a 
ramp  voltage,  the  value  of  which  is  displayed  on  a  strip  chart 
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A  current  probe  monitors  the  transducer  current  and  the 


output  is  displayed  on  an  oscilloscope. 

The  driving  amplitude  is  increased  until  a  cavitation  event  is 

detected.  Two  methods  of  detection  were  employed.  The  first  is  audio 

detection*  Vaporous  cavitation  events  are  almost  always  associated  with 

a  loud  audible  "click"  that  is  produced  by  a  collapsing  cavity.  Thus, 

cavitation  can  be  detected  by  the  operator  who  monitors  the  system 

through  a  headphone  set.  A  second  type  of  detection  scheme  relies  upon 

the  production  of  sonoluminescer.ee.  Sonoluminescence  is  an  emission  of 

light  associated  with  cavitation  events  In  liquids.  Its  origin  is  not 

clearly  understood,  but  recent  experimental  evidence  Indicates  that  it 

is  associated  with  the  deexcitation  and  radiative  recombination  of 
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chemical  species  formed  during  the  collapse  of  a  bubble.  mis 
emission  is  detected  with  a  photomultiplier  tube  mounted  underneath  the 
flask.  The  photomultiplier  signal  is  processed  by  a  noise  reduction 
system  to  reduce  the  occurrence  of  false  triggers.  When  the  detector 
senses  that  sonoluminescence  has  been  produced,  it  resets  the  ramp 
voltage  and  the  system  is  quiescent  for  a  predetermined  time  before  the 
ramping  phase  is  reinitiated.  During  the  quiescent  period  the  sound 
field  in  the  flask  is  generally  attenuated  by  60  dB  below  the  level 
associated  with  a  ramp  voltage  of  zero.  The  ramping/quiescent  cycle  is 
repeated  until  the  desired  number  of  data  points  are  taken. 

A-4  Experimental  Procedure 

The  discussion  to  follow  is  a  step-by-step  description  of  the 
procedure  used  to  acquire  a  typical  set  of  data.  An  effort  was  made  to 


follow  the  same  procedure  as  closely  as  possible  for  each  set  of  data. 
This  discussion  will  be  divided  into  two  sections:  one  for  the  sample 
preparation  and  the  other  for  data  acquisition. 

A-4.1  Procedure  for  Sample  Preparation 

The  water  used  for  these  measurements  was  acquired  from  a  bulk 
source  of  distilled  water  and  stored  in  a  25  liter  plastic  bottle.  In 
an  effort  to  eliminate  the  effect  of  any  variation  of  the  properties  of 
the  distilled  water  obtained  from  this  source  at  different  times,  the 
measurements  for  a  given  ionic  compound  were  all  made  with  water 
obtained  at  the  same  time.  The  water  used  for  measurement  of  the 
cavitation  threshold  at  a  given  concentration  was  all  drawn  from  the 
storage  bottle  at  the  same  time,  and  after  being  prefiltered  through  a 
teflon  filter,  was  stored  in  indi/idual  i  liter  teflon  PFA  bottles.  The 
pore  size  of  the  filter  was  chosen  to  be  the  same  as  that  of  the  in-line 
filter. 

The  Ionic  concentration  of  a  sample  was  adjusted  by  adding  a 
specified  volume  of  a  0.40  molar  stock  aqueous  solution  of  the  compound. 
The  stock  solution  was  made  by  adding  a  known  amount  of  the  compound 
(laboratory  grade)  to  one  liter  of  filtered  distilled  water  having  a 
specific  conductance  of  less  than  one  micromho/cm.  The  concentration  of 
the  stock  solution  was  calculated  from  the  measured  value  of  the 
specific  conductance.  The  conversion  was  made  according  to  the 
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where  is  the  ionic  concentration  in  moles/1,  A  is  the  molar 

2  -1 

conductance  in  mho  cm  mol  ,  and  K  is  the  specific  conductance  in 
mho/cm.  Small  amounts  of  the  compound  were  added  ,  if  necessary,  to 
bring  the  concentration  of  the  stock  solution  to  0.40  molar.  The  stock 
solution  was  kept  in  a  sealed  glass  flask  and  stored  in  a  closed  cabinet 
when  not  needed.  A  teflon  coated  magnetic  stirring  bar  was  kept  in  the 
solution  so  the  solution  could  be  mixed  immediately  before  being 
introduced  into  the  sample.  The  preparation  of  the  sample  will  now  be 
discussed. 

Assume  that  the  sample  preparation  system  contains  a  previously 
used  sample.  The  old  sample  is  pumped  out  of  the  system.  After  as  much 
of  the  liquid  as  possible  hag  been  pumped  out  of  the  system,  the  system 
is  pressurized  with  argon.  This  pressurization  forces  out  any  remaining 
liquid  except  for  the  portion  contained  within  the  quartz  flask  (about 
0.75  liters).  After  the  system  has  been  flushed,  about  1.75  liters  of 
prefiltered  distilled  water  is  introduced  into  the  system.  This  water 
is  circulated  for  about  15  minutes,  after  which  time  the  system  is 
flushed  once  again.  Then,  one  liter  of  prefiltered  distilled  water  is 
put  into  the  system.  This  water  along  with  that  remaining  in  the  flask 
is  the  water  used  for  the  next  data  run. 

Circulation  is  started  and  the  water  is  routed  through  the  organic 
and  ionic  removal  columns.  The  deionization  process  is  continued  for 
thirty  minutes  at  which  time  the  water  is  routed  so  as  to  bypass  the 
columns.  (However,  whenever  the  water  is  circulating,  it  passes  through 
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the  in-line  filter.)  Three  small  samples  of  the  deionized  water  (about 
5  ml  each)  are  removed  from  the  system.  The  conductance  of  these 
samples  is  measured  and  must  be  less  than  1.5  micromho/cm  in  order  for 
the  preparation  phase  to  progress  farther.  In  every  instance  the 
samples  always  met  this  requirement  after  thirty  minutes  of 
deionization.  In  most  cases  the  average  specific  conductance  was  less 
than  1.0  micromho/cm.  At  this  point  the  stock  solution  is  added.  In 
order  to  keep  the  system  volume  approximately  constant  for  all 
measurements,  an  additional  amount  of  water  equal  to  the  volume  of 
solution  to  be  added  was  removed  from  the  system.  The  solution  Is  then 
introduced.  Throughout  this  process  the  water  has  been  continuously 
circulated. 

After  the  solution  has  been  introduced,  the  degassing  phase  begins. 
The  absolute  pressure  during  degassing  was  maintained  at  0.30  bar,  and 
therefore,  the  dissolved  gas  pressure  of  argon  In  the  water  was  0.30 
bar.  As  described  in  reference  45,  the  system  has  been  calibrated 
against  actual  measurements  of  gas  content  made  with  a  Van  Slyke 
apparatus.  The  degassing/argon  saturation  phase  is  accompanied  by  a 
temperature  control  phase.  At  the  end  of  the  degassing  phase,  the 
temperature  of  the  sample  is  about  23.15°C.  Circulation  Is  stopped  and 
the  system  is  returned  to  atmospheric  pressure.  Three  more  samples  are 
drawn  for  the  purpose  of  determining  the  dissolved  ion  concentration  and 
surface  tension  of  the  sample.  The  system  is  allowed  to  rest  for  one 
hour.  This  rest  period  accomplishes  two  things.  It  allows  any  large 
bubbles  produced  by  the  circulation  and  agitation  to  dissolve,  or 
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stabilize,  or  rise  to  the  surface  of  the  liquid.  It  also  allows  the 
temperature  of  the  system  to  stabilize.  At  the  end  of  the  rest  period 
the  temperature  of  the  sample  in  the  quartz  flask  is  23.00  +  0.10°C. 
During  the  data  acquisition  phase,  the  temperature  of  the  sample  is 
maintained  at  23.00  +_  0.20°C.  If  the  temperature  of  the  sample  falls 
outside  this  range,  data  acquisition  is  stopped.  The  data  acquisition 
phase  will  be  described  next. 

A-A.2  Procedure  for  Data  Acquisition 

Near  the  end  of  the  one  hour  rest  phase,  the  acoustic  signal 
generation  system  is  calibrated  (described  in  detail  in  reference  A3). 
At  the  end  of  the  rest  period  the  data  acquisition  phase  begins.  The 
ramp  ic  initiated.  When  the  operator  audibly  detects  that  a  transient 
cavitation  event  has  occured ,  the  drive  speed  of  the  stripchart  recorder 
is  increased.  In  this  way  the  "audible"  cavitation  threshold  is  noted. 
If  the  audible  event  is  not  accompanied  by  soroluminescence ,  the  ramp 
continues  until  sonoluminescence  Is  detected.  When  it  is,  the  ramp 
voltage  is  reset,  the  chart  speed  is  reduced,  and  the  system  is 
quiescent  for  a  five  minute  period.  This  quiescent  phase  is  necessary 
to  allow  the  temperature  of  the  sample  to  stabilize  and  any  remnants  of 
the  cavitation  event  to  dissolve.  This  procedure  is  repeated  until 
either  15  data  points  are  taken  or  the  temperature  of  the  system  fails 
to  return  to  within  0.20°C  of  23.00°C  after  the  quiescent  period. 

The  stripchart  voltages  are  converted  to  pressure  amplitudes 
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according  to  the  following  relationship 
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PA  -  0.91  +  2.87  V  bat  (A2) 
A  s 

where  is  the  stripc.hart  voltage  corresponding  to  the  cavitation 
event.  (The  origin  of  this  relationship  is  described  in  reference  43.) 
Using  this  relationship  two  cavitation  thresholds  are  obtained:  an 
audible  threshold  and  a  sonoluminescence  threshold.  The  results  of 
these  measurements  are  presented  in  the  following  section. 

A-5  Results  of  the  Experiment 

Six  different  compounds  were  used  for  the  measurements:  NaF,  NaCl, 
Nal,  KCI,  KI,  and  KOH.  The  audible  threshold  as  a  function  of 
concentration  for  KI  and  Nal  were  presented  in  figures  18a  and  lflb.  The 
data  were  normalized  to  the  threshold  at  the  lowest  concentration  to 
facilitate  comparison  among  various  data  sets.  The  value  of  the 
threshold  at  the  lowest  concentration  is  typically  about  4.5  j-_  1  bar; 
such  variations  are,  unfortunately,  typical  in  cavitation  threshold 
measurements.  The  normalization  is  justifiable  since  it  is  the 
variation  in  the  threshold  that  is  sought  and  not  the  absolute  value. 
The  error  bars  in  figures  18a  and  18b  represent  plus-or-minus  one 
standard  deviation. 

The  lowest  value  of  the  concentration  was  determined  in  the 
following  manner.  The  sample  was  prepared  in  the  usual  way  except  that 
no  ionic  solution  was  introduced.  That  is,  the  threshold  of  the 
deionized  water  was  measured.  The  deionized  water  had  a  non-zero 
specific  conductance,  typically  about  1  micromho/cm.  In  order  to  assign 
a  concentration  to  this  deionized  sample,  it  was  assumed  that  the  non- 
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zero  conductance  was  due  to  the  presence  o£  whatever  compound  was  under 
investigation  at  that  time.  This  is  probably  not  a  valid  assumption, 
but  a  better  method  was  not  apparent.  Therefore,  the  values  of  the 
lowest  concentrations  are  to  he  considered  as  upper  limit’'  The  actual 
concentration  of  the  specific  compound  under  investigation  is  probably 
much  less  than  the  assigned  value. 

As  stated  above  our  apparatus  made  it  possible  to  measure  both  an 
audible  and  sonoluminescence  threshold.  In  figure  21a,  both  the  audible 
and  sonoluminescence  threshold  are  plotted  as  a  function  of  ion 
concentration.  In  every  Instance,  the  audible  threshold  is  lower  than 
the  sonoluminescence  threshold.  This  trend  was  present  in  all  six 
compounds  tested.  However,  the  variation  with  ion  concentration  was  the 
same  for  each  threshold  as  3hown  in  figure  21b.  That  is,  when  the 
thresholds  are  normalized  the  difference  between  them  is  negligible. 
This  trend  was  also  apparent  in  all  six  compounds.  From  this  dicussion 
it  may  be  concluded  that  for  measurements  of  the  variation  of  the 
transient  cavitation  threshold,  either  of  the  thresholds  are  applicable. 
Emphasis  has  been  placed  on  the  words  "variation"  and  "transient".  The 
absolute  threshold  does  not  appear  to  be  the  same  for  both  cases. 
Measurements  made  with  single  photon  counters  indicate  that  a  clear- 
cut  threshold  for  sonoluminescence  emission  from  stable  (gaseous) 
cavities  may  not  exist.  That  is,  pulsating  gas-filled  bubbles  may  emit 
light  even  when  the  pressure  amplitude  of  the  applied  acoustic  field  is 
extremely  small.  However,  when  transient  cavitation  occurs,  usually,  no 
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